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Semiconductor  manufacturing  is  one  of  today's  most  complicated  manufacturing 
processes.  Most  of  the  recent  research  in  semiconductor  manufacturing  focuses  on 
demand  satisfaction,  maximization  of  equipment  utilization,  and  maximization  of 
throughput  with  some  capacity  constraints.  In  this  research,  the  first  chapter  presents  an 
overview  of  the  semiconductor  manufacturing  technology  and  describes  the 
manufacturing  process.  It  also  discusses  a real  world  production  planning  application  in  a 
major  U.S.  semiconductor  company  and  proposes  a new  methodology  for  capacity 
planning. 

The  second  chapter  addresses  the  assignment  problem  of  wafers  to  alternative 
machine  groups  and  the  determination  of  equipment  requirements  by  recognizing 

viii 


capacity  limitations  of  the  individual  machines  as  well  as  reducing  operating  and 
investment  costs  related  to  the  machines.  A mathematical  model,  which  is  a variation  of 
the  well-known  generalized  assignment  problem  with  continuous/integer  variables  and 
fixed  charges,  is  developed  for  this  problem.  Given  the  intractability  of  the  model,  lower 
bounding  and  heuristic  methods  to  solve  the  problem  are  proposed. 

The  third  chapter  extends  the  problem  described  in  chapter  2 to  a dynamic 
environment  where  demands  vary  between  time  periods  and  inventory  carrying  is 
allowed.  We  investigate  the  tradeoff  between  allocating  a constant  capacity  for  the 
planning  horizon  and  holding  inventory  and  we  present  a mixed-integer  mathematical 
model  to  minimize  operating  and  investment  costs  related  to  the  machines  and  inventory 
holding  costs.  We  propose  a Lagrangian  relaxation  based  heuristic  approach  to  solve  the 
problem. 

The  fourth  chapter  presents  an  effective  algorithm  to  identify  a set  of  component 
families  in  order  to  effectively  to  sequence  individual  printed  circuit  boards  to  multiple, 
identical  insertion  machines.  The  objective  of  the  component  grouping  problem  is  to 
maximize  the  sum  of  similarities  between  each  component  type  and  a designated 
component  type  chosen  as  the  group  median  in  an  attempt  to  minimize  the  number  of 
machines  a printed  circuit  board  needs  to  visit.  Once  the  component  families  are  formed, 
two  objectives  are  considered  in  scheduling  the  printed  circuit  boards:  minimizing  the 
maximum  completion  time  (makespan)  and  the  mean  flow  time.  Finally,  in  the  fifth 
chapter  we  discuss  the  implications  and  conclusions  of  this  study  and  give  directions  for 
future  research. 
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CHAPTER  1 
INTRODUCTION 

1.1.  Semiconductor  Manufacturing 

This  chapter  presents  an  overview  of  the  semiconductor  manufacturing 
technology,  describes  the  manufacturing  process  introducing  the  terminology  that  will  be 
used  in  the  following  chapters.  It  also  discusses  a real  world  production  planning 
application  in  a major  U.S.  semiconductor  company  and  proposes  a new  methodology  for 
capacity  planning. 

1.1.1.  Semiconductor  Industry  Overview 

The  semiconductor  industry  is  an  excellent  example  of  growing  and  innovative 
industry.  Almost  non-existent  40  years  ago,  this  industry  has  grown  steadily  since  1985 
and  particularly  fast  since  1992,  with  annual  growth  exceeding  30  percent.  Industry  wide 
revenues  of  $27  billion  in  1984  and  $142  billion  in  1996  were  reported.  Projections  are 
that  these  will  be  $425  billion  by  2001.  This  rapid  growth  of  the  semiconductor  industry 
has  been  fueled  by  many  factors,  including: 

• A PC  market  that  continues  to  grow  - In  1995,  the  demand  for  PCs  resulted 
in  an  estimated  38  percent  of  chip  market  consumption. 
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• The  use  of  semiconductors  in  a wider  variety  of  products  - Semiconductors 
are  now  used  in  products  such  as  electronic  equipment,  appliances,  cellular 
phones,  and  automobiles.  Consumers  use  chips  embedded  in  goods  ranging 
from  the  microwave  oven  to  the  garage-door  opener. 

• The  increasing  microelectronics  content  of  electronics  products  - Chip 
content  of  the  average  electrical  product  was  4.3  percent  in  the  1970s,  6.8 
percent  in  the  1980s,  and  10.2  percent  in  the  early  1990s. 

• Emerging  international  markets,  including  Korea,  Taiwan,  and  China  - 
Today,  these  regions  represent  20  percent  of  world  consumption,  compared 
to  7 percent  a decade  ago. 

Such  a rapid  growth  in  the  semiconductor  industry  requires  a continuous 
expansion  of  the  manufacturing  infrastructure.  Much  of  the  current  restructuring  of  the 
U.S.  semiconductor  industry  may  be  understood  in  terms  of  efforts  to  reduce  cost  while 
producing  more  powerful  chips.  The  success  in  these  efforts  has  been  largely  based  upon 
innovation  and  technology  developments  such  as  the  implementation  of  more  advanced 
photolithography  techniques  and  the  utilization  of  larger  diameter  silicon  wafers. 

However,  the  productivity  gains  based  only  on  product  and  process  developments 
will  not  be  sufficient  to  offset  the  vigorous  growth  in  manufacturing  costs.  In  terms  of 
semiconductor-fabrication  equipment,  the  catch  up  between  new  manufacturing 
techniques  and  the  development  of  corresponding  equipment  has  clearly  become  a never- 
ending  race  in  which  each  surpasses  the  other.  Manufacturers  have  been  confronting  price 
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increases  of  28  percent  per  year  in  machinery.  This,  in  turn,  is  raising  the  high  cost  of 
building  new  state-of-the-art  semiconductor  facilities  and  investing  on  new  equipment. 

The  Semiconductor  Industry  Association  reports  that  new  fabrication  plants 
(called  “fabs”)  cost  an  average  of  $2.5  billion.  Within  the  next  two  years,  the  biggest  fabs 
will  cost  more  than  $5  billion.  Once  built,  fabs  can  become  obsolete  in  3 years  to  5 years. 
In  the  meantime,  market  demand  for  chips  can  change  rapidly,  making  it  difficult  for 
manufacturers  to  adjust  production  schedules  to  market  demands.  This  problem  is  being 
partially  alleviated  by  the  fact  that  chipmakers  are  now  tapped  into  networks  that  link 
them  to  their  PC  customers’  inventory  systems.  As  a result,  they  can  use  just-in-time 
manufacturing  techniques  to  match  supply  with  demand. 

However,  with  line  widths  (the  widths  of  individual  electronic  components  etched 
onto  the  wafer)  now  being  cut  to  fractions  of  a micron,  chipmakers  still  may  face 
diminishing  returns  on  their  next  generation  of  fab  manufacturing  investments.  Although 
the  cost  of  a fabrication  plant  doubled  between  1984  and  1990,  chipmakers  were  still  able 
to  double  the  number  of  transistors  on  a chip.  With  today’s  shrinking  line  widths, 
however,  it  will  be  difficult  for  the  next  generation  of  fabs  to  double  chip  performance 
again.  While  focusing  on  the  development  and  deployment  of  new  technologies, 
semiconductor  companies  must  simultaneously  embrace  other  methods  to  increase  their 
productivity  and  to  remain  globally  competitive.  These  methods  must  involve 
improvement  of  capital  utilization  through  the  increase  of  overall  equipment  effectiveness 
and  reduction  of  cycle  times  and  work-in-process  inventories. 
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1.1.2.  Semiconductor  Technology  Overview 

Semiconductors,  or  silicon  components,  are  the  basic  building  blocks  of  every 
piece  of  electronic  equipment.  Essential  in  computers,  they  are  also  part  of  everything 
from  microwave  ovens  to  automobiles.  The  components  themselves  are  small,  flat  chips 
of  silicon  with  miniaturized  electronic  circuits  etched  onto  them.  The  element  silicon  is 
the  most  commonly  used  semiconductor.  When  impurities  are  added  to  it  (a  process 
called  doping),  silicon  can  become  a conductor  or  insulator,  depending  on  the  dopant 
used.  The  term  semiconductor  is  also  used  in  a broader  sense  today  to  refer  to  electronic 
components  fabricated  from  semiconductor  materials,  such  as  integrated  circuits, 
processors,  and  auxiliary  chips. 

The  most  common  natural  semiconductors  used  in  electronics  are  the  elements  of 
silicon  and  germanium.  Silicon,  one  of  the  elements  of  common  beach  sand,  is  the  most 
frequently  used  semiconductor  material  because  it  can  be  used  at  a wide  variety  of 
temperatures.  Silicon  is  also  the  second  most  abundant  element  on  earth,  after  oxygen. 
Semiconductors  can  be  classified  according  to  their  material  (silicon,  germanium,  etc.), 
their  product  groups  (discrete  devices,  integrated  circuits,  etc.),  their  technology  (bipolar, 
MOS),  their  capability  to  handle  power,  and  their  customization  (standard,  custom  and 
semicustom).  We  will  be  focusing  on  the  classification  of  integrated  circuits  by  their 
technology. 

In  the  semiconductor  industry,  the  processes  to  make  integrated  circuit  (IC) 
products  may  vary  from  one  company  to  another,  yet  there  are  usually  some  strong 
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similarities.  By  design  technology,  there  are  only  a few  major  groups  or  divisions  such  as 
MOS  (Metal  Oxide  Semiconductor)  and  Bipolar.  MOS  transistors  are  small  and  consume 
relatively  little  power,  but  historically  were  relatively  slow.  Circuits  built  from  bipolar 
transistors  are  larger  and  faster,  but  consume  a great  deal  more  power.  Over  the  past 
decade,  MOS  technology  became  more  sophisticated,  resulting  in  increased  processor 
speed. 

These  two  divisions  have  usually  subdivisions.  CMOS  (Complementary  MOS)  is 
now  the  most  frequently  used  technology  in  systems  ranging  from  PCs  to  mainframes  to 
supercomputers.  Each  of  these  subdivisions  has  its  own  method  for  making  its  products. 
Within  each  subdivision,  many  products  will  be  manufactured.  An  ECL  (Emitter  Coupled 
Logic)  group  will  not  produce  just  one  type  of  product,  in  fact,  they  will  possibly  make 
hundreds  of  different  products. 

In  semiconductor  manufacturing,  a combination  of  numbers  and  letters  is  used  to 
identify  a particular  product  type  that  is  being  made.  Similar  to  the  automobile  industry 
where  companies  (e.g.  Chrysler)  have  different  divisions  (e.g.  Dodge,  Plymouth,  Jeep) 
that  manufacture  cars  with  different  names  (e.g.  Stratus,  Neon,  Cherokee),  computer 
chips  in  semiconductor  industry  have  number/letter  codes  as  illustrated  in  Figure  1.1.  The 
major  difference  between  the  automobile  industry  and  semiconductor  industry  is  that 
General  Motors  does  not  manufacture  Neons  nor  does  Chrysler  manufacture  Cavaliers. 

On  the  other  hand,  in  semiconductor  industry  many  companies  may  be  making  the  same 


XYZ90  product. 
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Bipolar  MOS 


ECL  TTL  CMOS 

XYZ90  123AB  7SX  Logic 


ECL:  Emitter  Coupled  Logic 
TTL:  Transistor-Transistor  Logic 
CMOS:  Complimentary  Metal  Oxide  Semiconductor 
HMOS:  High  density  Metal  Oxide  Semiconductor 
NMOS:  N-channel  Metal  Oxide  Semiconductor 


Figure  1.1.  Classification  of  integrated  circuits1 


HMOS  NMOS 

5Y  Mem  N2-7X1 


1.1.3.  Semiconductor  Manufacturing  Process 

Manufacturing  a computer  chip  is  a complex  process  involving  hundreds  of  steps 
and  requiring  from  a few  days  to  up  to  three  months  of  processing  time  to  complete. 
Chips  require  incredible  precision,  careful  attention  to  detail,  and  constant  testing.  The 
manufacturing  process  can  be  divided  into  five  stages:  the  mask  design,  the  crystal 
growth  and  slice  preparation,  the  fabrication  of  the  device,  the  assembly  and  the  final 
testing. 


1 All  integrated  circuit  names  are  fictitious. 
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OXIDATION 


DIFFUSION 


DEPOSITION  MASKING 


Initial  Oxidation  N+  Diffusion 


1st  Deposition  1st  Mask 


Boron  Oxidation  Isolation  Diffusion  2nd  Deposition  2nd  Mask 


3rd  Oxidation 


P+  Diffusion 


3rd  Deposition  Metal  Mask 


Figure  1.2.  Semiconductor  manufacturing  operation  steps 


In  IC  manufacturing,  each  letter  and  number  code  represents  a distinct  product 
and  each  product  has  its  own  special  set  of  operations  such  as  “Oxidation,”  “Diffusion,” 
“Masking,”  etc.  which  must  be  performed  to  make  it  successfully  function.  Some  of  the 
work  performed  in  these  operations  are  shown  in  Figure  1.2.  These  processing  steps  are 
combined  into  a certain  work  flow  order  when  the  wafers  are  processed  and  each 
particular  sequence  throughout  its  fabrication  process  visiting  different  work  centers. 
Specifically,  each  wafer  makes  multiple  visits  to  the  same  work  center  at  different  points 
in  the  fabrication  process.  The  common  practice  in  the  layout  of  facilities  for  different 
types  of  processes  is  to  group  similar  operations  together:  furnaces  are  grouped  in  one 
area,  ion  implanters  in  another,  and  so  on.  This  layout  requires  wafers  moving  back  and 
forth  between  work  centers  but  it  allows  the  utilization  of  the  same  equipment  to  process 
wafers  at  different  operation  steps.  For  instance,  all  masking  operations  may  be 
completed  on  the  same  machine  or  in  the  same  work  center  although  the  wafers  are 
required  to  travel  to  other  processing  areas  between  these  masking  operations. 
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Table  1.1  illustrates  an  example  of  operation  naming  and  sequencing  practice 
along  with  the  production  volume  of  each  product  type2.  Note  that  the  operation  steps 
may  be  same  or  distinct  for  each  product  type  and  the  operation  name  indicates  both  the 
product  type  and  the  process  to  be  performed  on  that  product  type.  The  product  and 
operation  incidences  with  the  scheduled  volume  of  products  at  each  operation  step  are 
also  given  in  Table  1.2  as  an  example.  This  example  depicts  the  case  of  three  types  of 
wafer  and  five  types  of  process  they  require  during  the  manufacturing  process:  oxidation, 
masking,  deposition,  diffusion,  and  ion  implantation.  These  processes  are  repeated  at 
different  steps  and  there  are  21  distinct  operation  steps.  The  operation  name  describes  the 
type  of  the  process  at  each  step  and  the  wafer  type  to  be  processed.  For  example,  {CM- 
0X1 -X}  can  be  interpreted  as  the  first  oxidation  step  on  a X type  CMOS  wafer.  The 
number  of  CMOS/X  wafers  to  be  processed  each  day  at  this  stage  is  80.  Similarly,  {CM- 
DEP1-XYZ}  denotes  the  first  deposition  process  for  all  X,  Y,  and  Z types  of  CMOS 
wafers.  The  scheduled  volume  at  this  operation  step  is  300  wafers,  which  is  the  sum  of 
the  volumes  of  CMOS/X,  CMOS/Y,  and  CMOS/Z  chips,  since  they  all  require  the  same 
deposition  process  and  can  be  processed  with  the  same  setup.  Assuming  100  % yield  per 
operation  in  this  example,  the  number  of  wafers  per  day  remain  constant  through  the 
manufacturing  process. 


2 The  names  are  for  illustrative  purposes  and  fictitious. 
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Table  1.1.  Operation  names  and  sequences  by  products 


PRODUCT 

VOLUME 

OPERATION 

OPERATION 

NAME 

(Wafers/Day) 

NUMBER 

NAME 

CMOS  /X 

80 

1 

CM-OX1-X 

4 

CM-M1-XY 

6 

CM-DEP1-XYZ 

7 

CM-DIFF1-X 

10 

CM-OX2-X 

12 

CM-M2-XZ 

14 

CM-IMP1-XZ 

16 

CM-DEP2-XY 

18 

CM-DIFF2-XY 

20 

CM-IMP2-XZ 

CMOS  / Y 

100 

2 

CM-OX1-Y 

4 

CM-M1-XY 

6 

CM-DEP1-XYZ 

8 

CM-DIFF1-Y 

11 

CM-OX2-YZ 

13 

CM-M2-Y 

15 

CM-IMP1-Y 

16 

CM-DEP2-XY 

18 

CM-DIFF2-XY 

21 

CM-OX3-Y 

CMOS  /Z 

120 

3 

CM-OX1-Z 

5 

CM-M1-Z 

6 

CM-DEP1-XYZ 

9 

CM-DIFF1-Z 

11 

CM-OX2-YZ 

12 

CM-M2-XZ 

14 

CM-IMP1-XZ 

17 

CM-DEP2-Z 

19 

CM-MM-Z 

20 

CM-IMP2-XZ 

Table  1.2.  Product  / Operation  incidences  and  scheduled  volumes 
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1 .2.  A Production  Planning  Application 

This  section  presents  a production  planning  application  in  a major  U.S. 
semiconductor  company  that  develops,  manufactures,  and  markets  semiconductor 
technologies,  products,  packaging  and  services  for  their  internal  and  external  customers. 
The  manufacturing  facilities  range  from  low  volume,  high  product  variation  research  and 
development  lines  to  high  volume  lines  with  both  low  and  high  volume  variation. 

The  manufacturing  facility  where  this  study  has  been  conducted  has  two  wafer 
fabrication  areas.  These  fabs  have  been  connected  to  each  other  recently  to  facilitate  the 
routing  of  work-in-process  (WIP)  inventories  and  provide  possible  use  of  additional 
machine  capacity.  The  total  number  of  machines  in  these  fabs  is  around  1200.  The  cost  of 
a machine  varies  roughly  between  $2  million  and  $7  million  and  the  facility  purchases  50 
machines  on  the  average  each  year  to  replace  the  obsolete  machines  or  to  add  capacity  to 
meet  demand  forecasts.  Hence,  capacity  planning  is  a strategic  decision  since  the  cost  of 
investing  in  new  equipment  is  significantly  high. 

1.2.1.  Machine  Planning  System  Overview 

The  Machine  Planning  System  (MAPS)  is  designed  to  facilitate  manufacturing 
capacity  planning,  which  relates  to  estimating  the  amount  and  type  of  process  equipment 
required  to  produce  a given  product  mix  over  a particular  time,  and  using  these  estimates 
for  current  managerial  decisions  and  new  facility  construction.  When  only  a unique 
machine  can  complete  a certain  operation  this  calculation  is  relatively  straight  forward; 
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simply  sum  the  number  of  operations  to  be  processed,  divide  by  the  production  capacity 
of  the  machine  type,  its  yield,  and  its  utilization  to  arrive  at  the  number  of  machines 
required.  However,  in  most  manufacturing  facilities  there  are  often  several  different 
machine  sets  that  can  be  used  to  process  a particular  operation.  In  the  case  where  no 
single  machine  group  has  the  sufficient  capacity  to  process  all  the  required  volume,  the 
production  WIP  needs  to  be  allocated  among  alternate  machine  groups  and  the  capacity 
across  all  possible  machine  groups  is  assessed. 

Each  facility  in  this  company  uses  its  own  internally  developed  Machine  Planning 
System  and  each  system  is  linked  to  others’  databases  (However,  the  fabrication  areas 
remain  independent).  The  major  drawbacks  of  the  present  system  are  the  incompatibility 


Operations 


Primary  Machine  Group 


CM-OX1-X 

CM-OX1-Z 


Operations 


Secondary  Machine  Group 


CM-OX1-Y 

CM-OX2-YZ 


Operations 


Tertiary  Machine  Group 


CM-OX2-X 

CM-OX3-Y 


Figure  1.3.  Cascades  in  MAPS 
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of  one  system  with  another,  which  causes  difficulties  in  data  integration  and  exchange, 
and  the  age  which  results  in  extensive  and  expensive  maintenance  and  upgrade 
requirements.  The  machine  planning  engine  in  MAPS  deals  with  groupings  of  similar 
machines  called  “machine  groups”.  Special  names  are  assigned  to  different  machine 
groups  and  all  machines  within  a group  are  equivalent.  Machine  groups  are  allowed  to  be 
dedicated  to  different  operations  but  not  individual  machines  within  a group. 

1.2.2.  Cascading  in  MAPS 

Cascading  is  the  term  used  for  allocating  production  work  across  alternate  parallel 
machine  groups  in  capacity  planning.  Given  a manufacturing  line,  there  is  typically  a mix 
of  multiple  machine  types  with  different  production  characteristics  that  can  perform  the 
same  operation.  One  machine  type  is  generally  preferred  over  another,  though  this 
preference  may  differ  according  to  either  the  process  or  product  involved.  As  a result, 
priorities  are  assigned  to  each  unique  process/product  which  we  define  as  an  operation 
and  work  is  allocated  across  alternate  machine  groups  on  that  basis. 

The  machine  groups  in  MAPS  are  classified  as  primary  machine  group,  secondary 
machine  group,  tertiary  machine  group,  and  so  on,  for  any  operation  set  where  primary 
machines  are  the  most  efficient  machines  that  can  process  the  specified  operations  and 
secondary,  tertiary,  etc.  machines  are  the  alternative  machines  in  case  additional  machine 
capacity  is  needed.  Cascading  in  MAPS  refers  to  the  allocation  of  secondary  machine 
capacity  to  operations  whose  primary  machine  capacity  has  been  exceeded. 
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In  the  current  system,  WIP  allocation  is  carried  out  sequentially.  In  general  terms, 
the  machine  groups  that  can  process  some  of  the  same  operations  (product/process 
combinations)  are  assigned  in  ordered  sequence  to  a cascade  group.  The  production 
volumes  for  all  operations  assigned  to  the  cascade  group  are  first  allocated  to  the  primary 
machine  group,  once  its  capacity  is  reached  the  remaining  volumes  are  allocated  to  the 
secondary  machine  group,  then  to  the  tertiary,  etc.,  until  all  machine  groups  have  been 
allocated  with  as  much  work  as  their  defined  capacity.  If  unsolutioned  volumes  remain 
after  loading  all  machine  groups  in  the  cascade  groups,  the  remaining  work  is  allocated 
first  to  the  primary  machine  group,  for  all  operations  defined  on  that  machine  group,  if 
not  possible,  then  sequentially  to  the  secondary  and  tertiary  machine  groups,  etc.,  until  all 
production  volumes  have  been  assigned.  Thus,  if  the  primary  machine  group  is  defined 
with  the  capability  of  processing  all  operations  within  the  cascade  group,  all  unsolutioned 
volume  would  be  allocated  to  just  that  single  machine  group,  leaving  it  the  only  machine 
group  needing  additional  machines  to  contain  the  capacity  requirement.  The  term 
cascading  comes  from  the  analogy  of  describing  each  machine  group  as  a pool  of 
capacity,  the  production  volumes  flow  first  into  the  primary  pool,  once  it  is  full  the 
remaining  volume  (if  any)  flows  to  secondary  pool,  and  so  on,  like  a cascading  waterfall. 
A step  by  step  description  of  the  algorithm  is  as  follows: 

Current  Cascading  Algorithm: 

Step  0 : Compute  the  number  of  machines  required  for  each  operation. 
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Step  1 : Assign  all  unique  operations  to  machine  groups  that  are  capable  of 

performing  them. 

Step  2 : Fill  up  primary  machine  group  in  the  given  order  of  the  operations. 

Step  3 : Repeat  with  secondary  machines  if  primary  machines’  capacity  has  been 

exceeded. 

Step  4 : Repeat  until  all  cascaded  machines  are  used  up  or  all  operations  are 

assigned. 

Step  5 : Add  machines  to  satisfy  demand  (if  needed). 

Step  6 : Repeat  for  all  periods.. 

The  data  required  for  the  algorithm  include  daily  production  capacity,  yield,  and 
utilization  of  each  machine  type,  number  of  machines  of  each  type  in  each  machine  group 
and  scheduled  volume  per  day  of  each  product  for  the  specified  operation  steps.  The 
production  capacity  is  the  number  of  wafers  processed  per  day  and  is  implicitly  a function 
of  the  raw  processing  time,  setup  time,  and  operator  intervention  time.  The  yield  is  the 
ratio  of  number  of  good  wafers  at  the  end  of  the  process  to  the  number  of  wafers  that  go 
in.  Utilization  is  the  net  production  time  of  the  machine  group  during  the  24-hour  period. 
Using  this  data,  the  number  of  machines  required  by  each  operation  on  each  machine 
group  is  computed  as  follows  based  on  the  assumption  that  all  wafers  are  processed  on 
that  machine  group: 

Volume 


Number  of  machines  needed  = 


Production  Capacity  * Yield  * Utilization 
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1.2.3.  A New  Methodology  for  Cascading 

In  this  section  we  introduce  a new  cascading  methodology  that  includes 
enhancements  designed  to  improve  some  inherent  limitations  within  the  current  MAPS 
cascading  methodology.  Primary  among  these  is  the  restriction  that  all  operations  defined 
within  the  cascade  group  have  the  same  ordinal  allocation  priority  across  the  machine 
groups  in  a cascade  set.  Although  one  machine  group  is  often  preferred  over  another  for  a 
particular  operation,  this  machine  preference  may  vary  among  the  designated  operations 
due  to  some  factors  such  as  speed,  yield,  logistics,  robustness,  etc.  Thus,  a sequential 
WIP  flow  such  as  current  cascading  can  result  in  an  impractical  and  unrealistic  capacity 
consumption  solution.  The  new  cascading  algorithm  introduces  a declared  allocation 
priority  for  each  operation  across  the  machine  groups  within  the  cascade  set.  It  is  a more 
flexible  bi-directional  WIP  allocation  methodology  that  should  result  in  a more  practical 
machine  consumption  solution. 

At  the  time  that  an  operation  is  associated  with  a machine  group,  whether 
cascaded  or  not,  it  is  assigned  a load  priority  sequence  number.  This  priority  number  is  a 
non-unique  number  describing  a weighted  preference  for  the  order  in  which  volumes 
should  be  processed  through  the  machine.  A priority  value  of  1 is  the  highest  priority  or 
first  operation  to  be  processed,  a 9 would  be  the  lowest  priority  or  last  operation 
processed  on  that  machine  group3.  All  operations  with  the  same  priority  number  are 
loaded  in  the  original  sequence. 


3 A single  digit  is  used  for  illustrative  purposes. 
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Once  the  cascade  group  is  defined  with  sequentially  ordered  machine  groups,  an 
allocation  priority  number  is  assigned  defining  an  operation’s  preference  for  each  capable 
machine  group.  This  priority  is  a non-unique  number  describing  a weighted  preference 
for  where  operation  should  be  processed.  Factors  to  consider  include  speed,  yield, 
logistics,  robustness,  etc.  A priority  value  of  1 designates  the  most  preferred  machine 
group,  a 9 the  least  preferred  for  each  operation.4  If  multiple  machine  groups  are 
designated  the  same  preference  number,  operations  are  assigned  based  upon  the  original 
sequence  of  machine  groups  within  the  cascade  set. 

After  both  priorities  are  defined,  a composite  priority  is  created  for  every 
operation  on  each  machine  group  within  the  cascade  set  called  as  “cascade  priority.”  This 
is  a two-digit  integer  value  where  the  right  digit  is  the  first  priority  number  described 
above  and  the  left  digit  is  the  other.  If  the  operation  cannot  be  processed  by  a machine 
group,  the  restricted  priority  of  99  is  assigned. 

Based  on  these  cascade  priorities  the  proposed  methodology  assigns  operations 
across  machine  groups  in  two  stages.  In  the  first  stage,  the  operations  are  sorted  lowest  to 
highest  in  order  of  their  cascade  priority  number  for  each  machine  group.  At  the  first 
pass,  operations  whose  priorities  start  with  1 are  assigned  in  the  order  of  the  second  digit 
of  their  priorities  starting  with  the  first  machine  group  in  the  order  (primary  machine 
group  in  the  current  algorithm)  and  followed  by  the  second  machine  group  (secondary 
machine  group  in  the  current  algorithm),  then  third  machine  group  (tertiary  machine 


4 Again,  a single  digit  priority  is  used  for  illustrative  purposes. 
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group),  and  so  on.  Volumes  are  assigned  until  the  individual  machine  group  is  at  capacity 
for  the  number  of  available  machines.  Left  over  volumes  are  assigned  either  in  the 
following  passes  on  the  other  machine  groups  or  at  the  end  of  the  first  stage  as 
unsolutioned  operations. 

After  all  machine  groups  have  been  assigned  with  operations  whose  priorities  start 
with  1 , a second  pass  is  performed  in  the  same  manner  where  operations  whose  priorities 
start  with  2 are  assigned  in  the  order  of  the  second  digit  of  their  priorities.  This  allocation 
procedure  is  continued  until  all  volumes  are  assigned  or  all  machine  groups  in  the  cascade 
set  are  at  capacity  for  the  number  of  available  machines.  Any  remaining  unsolutioned 
operation  volumes  are  then  assigned  to  their  most  preferred  machine  group  (i.e.  to 
machine  group  with  the  lowest  cascade  priority  number). 

By  allocating  the  volumes  based  upon  the  cascade  priorities  it  is  possible  that  one 
or  more  machine  groups  in  the  cascade  set  be  loaded  over  capacity  while  other  machine 
groups  are  underloaded.  Since  this  situation  contradicts  the  primary  goal  of  minimizing 
additional  machine  investment  a second  stage  in  the  allocation  is  undertaken  for 
operations  on  overloaded  machine  groups  only.  Operation  volumes  are  reassigned  from 
overloaded  machine  groups,  based  on  reverse  cascade  priority,  to  underloaded  machine 
groups. 

In  what  follows  is  a step  by  step  description  of  the  proposed  cascading 
methodology.  A brief  example  is  presented  in  the  next  section. 

New  Cascading  Algorithm: 

Step  0 : Compute  the  number  of  machines  required  for  each  operation. 
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Step  1 : 


Step  2 : 


Step  3 : 


Step  4 : 


Step  5 : 


Step  6 : 


Step  7 : 


Assign  the  operations  to  the  first  machine  group  in  the  non-decreasing 
order  of  cascade  priorities  that  start  with  1,  until  all  machines  are  used  up 
or  all  operation  volumes  are  assigned. 

Repeat  with  the  second  machine  group  and  continue  until  the  last  machine 
group  is  considered. 

Return  to  the  first  machine  group  for  the  second  pass  (considering  cascade 
priorities  between  21-29  this  time)  and  repeat  the  same  procedure  with 
unassigned,  partially  or  completely,  operations.  Stop  after  the  procedure  is 
executed  for  a predetermined  number  of  passes.5 
Compute  machine  shortages  (if  any)  for  unsolutioned  operations  and 
assign  them  to  their  most  preferred  machine  groups. 

If  one  or  more  machine  groups  are  overloaded  while  there  exist  excess 
capacity  on  another  (other)  machine  group(s),  transfer  operation(s)  in  the 
non-increasing  order  of  the  cascade  priorities  scanning  all  underloaded 
machine  groups  in  reverse  cascade  priority. 

Continue  transferring  operations  until  overloaded  machine  group(s)  is 
(are)  reduced  to  full  capacity  or  underloaded  machine  group(s)  is  (are) 
loaded  up  to  full  capacity. 

Repeat  for  all  periods. 


5 The  number  of  passes  is  set  equal  to  the  number  of  machine  groups. 
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Step  8 : If  additional  machines  are  needed,  add  them  and  redo  assignments 

following  Step  1 through  Step  3. 

1.2.4.  Cascading  example 

This  example  describes  the  production  planning  of  3 machine  groups  and  6 
operations  for  2 periods  using  the  proposed  cascading  algorithm.  There  are  5 machines  in 
the  machine  group  1 and  3 in  machine  groups  2 and  3.  The  scheduled  processing  volumes 
of  wafers  per  day  are  400  for  all  operations  during  the  first  period  and  520,  500,  400,  400, 
500,  520,  respectively,  during  the  second  period.  The  yields  of  machine  groups  are  1.00, 
0.99,  0.98,  their  production  capacities  are  300,  325,  333,  and  their  utilizations  are,  0.90, 
0.80,  0.75,  respectively.  The  cascade  priorities  are  as  shown  on  Table  1.3.  The  first  digit 
is  the  preference  of  the  machine  group  across  all  machine  groups  for  the  operation  and  the 
second  digit  is  the  assignment  preference  of  the  operation  among  all  operations.  For 


Table  1.3.  Priority  table  for  the  cascading  example 


Oper  1 

Oper  2 

Oper  3 

Oper  4 

Oper  5 

Oper  6 

Machine  Group  1 

11 

12 

23 

99 

34 

35 

Machine  Group  2 

23 

24 

99 

11 

99 

32 

Machine  Group  3 

99 

23 

34 

35 

11 

12 
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Table  1.4.  Assignment  of  operations  at  each  pass 


Machine  Group  1 

Machine  Group  2 

Machine  Group  3 

1st  Pass  - 

— ► Operation  1 

Operation  2 

Operation  4 

Operation  5 
Operation  6 

2 Pass 

— ► Operation  3 

Operation  1 
Operation  2 

Operation  2 

3rd  Pass  - 

— ► Operation  5 

Operation  6 

Operation  6 

Operation  3 
Operation  4 

instance,  priority  of  operation  1 for  machine  group  1 is  1 1 , which  means  that  at  the  first 
pass  this  operation  will  be  considered  first  to  be  assigned  to  machine  group  1,  then 
operation  2,  which  has  priority  12,  will  be  considered  for  the  same  machine  group.  At  the 
same  first  pass,  operation  4 will  be  assigned  to  machine  group  2 (it  has  cascade  priority 
1 1)  and  operation  5 to  machine  group  3 (it  has  cascade  priority  1 1).  At  the  second  pass, 
operation  3 will  be  assigned  to  machine  group  1 if  that  machine  still  has  capacity, 
operations  1 and  2 will  be  assigned  to  machine  group  2 if  they  have  not  been  completely 
assigned  to  machine  group  1 and  any  remaining  WIP  of  operation  2 will  be  assigned  to 
machine  group  3.  The  same  procedure  will  be  applied  for  a third  pass.  Note  that  we  allow 
partial  assignment  of  an  operation  to  more  than  one  machine  group  depending  on  the 
machine  capacities.  The  assignment  order  in  this  procedure  at  each  pass  is  given  on  Table 


1.4. 
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Table  1.5.  Production  plan  after  the  first  stage 


Machine  Group  1 Machine  Group  2 Machine  Group  3 


Period  1 

Period  2 

Period  1 

Period  2 

Period  1 

Period  2 

Operation  1 

1.48 

1.93 

0.00 

0.00 

0.00 

0.00 

Operation  2 

1.48 

1.85 

0.00 

0.00 

0.00 

0.00 

Operation  3 

1.48 

1.48 

0.00 

0.00 

0.00 

0.00 

Operation  4 

0.00 

0.00 

1.55 

1.55 

0.00 

0.00 

Operation  5 

0.00 

0.00 

0.00 

0.00 

1.63 

2.04 

Operation  6 

0.24 

0.00 

0.00 

1.11 

1.37 

0.96 

Available 

5 

5 

3 

3 

3 

3 

Required 

4.68 

5.26 

1.55 

2.66 

3.00 

3.00 

Whole 

5 

6 

2 

3 

3 

3 

Delta 

0 

-1 

1 

0 

0 

0 

SHORTAGES 

Operation  3 

0.00 

0.26 

0.00 

0.00 

0.00 

0.00 

Table  1.5  illustrates  the  production  plan  and  additional  machine  requirements  at 
the  end  of  the  first  stage  after  three  passes  are  completed.  Notice  that  following  its 
cascade  priority  unsolutioned  WIP  of  operation  3 is  assigned  to  machine  group  1 , 
resulting  this  machine  group  to  be  overloaded  in  the  second  period.  Notice  also  that 
although  machine  group  2 is  under  capacity  it  cannot  process  operation  3 since  the 
cascade  priority  of  operation  3 for  machine  group  2 is  99. 
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Table  1.6.  Final  production  plan 


Machine  Group  1 

Machine  Group  2 

Machine  Group  3 

Period  1 

Period  2 

Period  1 

Period  2 

Period  1 

Period  2 

Operation  1 

1.48 

1.93 

0.00 

0.00 

0.00 

0.00 

Operation  2 

1.48 

1.59 

0.00 

0.28 

0.00 

0.00 

Operation  3 

1.48 

1.48 

0.00 

0.00 

0.00 

0.00 

Operation  4 

0.00 

0.00 

1.55 

1.55 

0.00 

0.00 

Operation  5 

0.00 

0.00 

0.00 

0.00 

1.63 

2.04 

Operation  6 

0.24 

0.00 

0.00 

1.11 

1.37 

0.96 

Available 

5 

5 

3 

3 

3 

3 

Required 

4.68 

5.00 

1.55 

2.94 

3.00 

3.00 

Whole 

5 

5 

2 

3 

3 

3 

Delta 

0 

0 

1 

0 

0 

0 

SHORTAGES 


The  final  production  plan  following  the  second  stage  is  shown  on  Table  1.6.  Unassigned 
WIP  of  operation  3 remains  assigned  to  machine  group  1,  its  most  preferred  machine 
group,  while  a portion  of  operation  2 is  transferred  from  machine  group  1 to  machine 
group  2 since  the  latter  has  unused  capacity.  Notice  that  transferred  amount  of  “0.26”  in 
machine  group  1 corresponds  to  “0.28”  in  machine  group  2 due  to  production  capacity 


difference  between  the  two  machine  groups. 
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1 .3.  The  Organization  of  this  Dissertation 

This  rest  of  this  dissertation  consists  of  three  chapters  with  each  chapter  being  an 
independent  essay  and  complete  within  itself.  The  first  two  chapters  address  the  problem 
of  operation/machine  type  assignment  and  machine  type  duplication  in  semiconductor 
manufacturing  as  an  introduction  is  presented  in  this  chapter.  Chapter  2 introduces  a 
mixed  integer  programming  model  of  the  problem  considering  the  operating  and 
procurement  costs  of  the  machines  and  provides  efficient  lower  bounding  and  heuristic 
methods  to  solve  the  problem.  Chapter  3 is  an  extension  of  this  problem  to  a multi-period 
production  planning  perspective  where  inventories  may  be  carried  over  periods.  Chapter 
4 proposes  an  effective  algorithm  to  identify  the  component  families  from  a large  number 
of  component  types  in  the  manufacture  of  printed  circuit  board  assembly  system  and  to 
schedule  the  boards  to  insertion  machines.  Finally,  Chapter  5 summarizes  this  study  and 


gives  directions  for  future  research. 


CHAPTER  2 

OPERATION/MACHINE  TYPE  ASSIGNMENT  IN  SEMICONDUCTOR 

MANUFACTURING 

2.1.  Introduction 

Manufacturing  a computer  chip  is  a complex  process  involving  hundreds  of  steps 
and  requiring  from  a few  days  to  up  to  three  months  of  processing  time  to  complete.  The 
manufacturing  process  consists  of  various  operation  types  such  as  oxidation,  deposition, 
lithography,  diffusion,  etching,  ion  implantation,  etc.  Some  of  these  operation  types  are 
repeated  several  times  to  build  different  layers  on  the  wafer.  Each  wafer  type  follows  a 
particular  sequence  throughout  its  fabrication  process  visiting  different  work  centers. 
Specifically,  each  wafer  makes  multiple  visits  to  the  same  work  center  at  different  points 
in  the  fabrication  process.  The  common  practice  in  the  layout  of  work  centers  for 
different  types  of  processes  is  to  group  similar  operation  types  together:  furnaces  are 
grouped  in  one  area,  ion  implanters  in  another,  and  so  forth.  This  layout  requires  wafers 
moving  back  and  forth  between  work  centers  but  it  allows  the  utilization  of  the  same 
equipment  to  process  wafers  at  different  steps.  For  instance,  all  oxidation  processes  may 
be  completed  on  the  same  machine  or  in  the  same  work  center  although  the  wafers  are 
required  to  travel  to  other  work  centers  between  these  oxidation  process  steps. 


25 


26 


The  capacity  allocation  and  machine  duplication  in  semiconductor  manufacturing 
involves  the  assignment  of  a set  of  wafers  requiring  certain  process  types  to  different 
machine  groups,  and  the  determination  of  equipment  requirements.  We  represent  the 
problem  of  assigning  operations  to  different  machine  groups  where  machine  duplication 
is  permitted  as  a variation  of  the  well  known  Generalized  Assignment  Problem  (GAP). 
Given  a set  of  tasks  to  be  assigned  to  a group  of  agents  and  the  cost  of  performing  each 
task  by  each  agent,  the  GAP  consists  of  minimizing  the  cost  of  assigning  tasks  to  agents 
such  that  each  task  is  assigned  exactly  to  one  agent  subject  to  the  availability  of  a single 
resource  type  that  is  consumed  by  the  agents  when  performing  these  tasks  (Ross  and 
Soland,  1975). 

The  applications  of  GAP  include  computer  job  assignments  in  computer 
networks,  assignment  of  software  development  tasks  to  programmers,  scheduling  variable 
length  commercials  into  time  slots,  scheduling  payment  on  accounts  where  lump  sum 
payments  are  prespecified  (Balachandran,  1972),  assignment  of  ships  to  yards  for 
overhaul  (Gross  and  Pinkus,  1972),  communication  network  design  models  with  certain 
node  capacity  constraints  (Grigoriadis,  Tang,  and  Woo,  1974).  Ross  and  Soland  (1977) 
show  how  to  formulate  and  solve  some  special  facility  location  problems  as  GAPs.  Fisher 
and  Jaikumar  (1981)  address  the  vehicle  routing  problem  based  on  a GAP  model  where 
the  vehicle  fleet  delivers  products  from  a central  warehouse  to  scattered  customers.  An 
extension  to  GAP  is  known  as  the  Multi-Resource  Generalized  Assignment  Problem 
(MRGAP)  where  the  agents  consume  multiple  resources  when  performing  their  tasks 
(Gavish  and  Pirkul,  1991).  MRGAP  was  used  in  the  design  of  distributed  computer 
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systems  (Gavish  and  Pirkul,  1986)  and  in  trucking  industry  (Murphy,  1986).  The 
interested  reader  is  referred  to  Cattrysse  and  Van  Wassenhove  (1992)  for  a broad  survey 
of  the  algorithms  developed  for  GAP  and  MRGAP. 

Both  GAP  and  MRGAP  assume  that  each  task  is  assigned  to  a single  agent. 
However,  in  semiconductor  industry,  as  well  as  in  many  other  industries,  it  is  a common 
production  planning  practice  to  assign  a task  to  two  or  more  agents.  In  this  chapter,  we 
discuss  an  operation/machine  type  assignment  problem  where  task  splitting  is  allowed 
and  resource  expansion  is  possible.  The  remainder  of  the  chapter  is  organized  as  follows. 
In  Section  2.2  we  describe  the  problem  and  formulate  a mathematical  model.  In  Section 
2.3  we  discuss  Lagrangian  relaxation  and  Lagrangian  decomposition  techniques  to  obtain 
lower  bounds  to  the  problem,  and  a solution  procedure  for  the  problem  is  described  in 
Section  2.4.  In  Section  2.5,  we  describe  experimental  analysis  of  the  proposed  solution 
method  on  both  randomly  generated  and  industry  data. 

2.2.  Problem  Description  and  Formulation 

We  model  the  problem  of  assigning  individual  operations  to  predetermined 
machine  groups  where  machine  duplication  is  allowed  as  a variation  of  GAP  with  both 
continuous  and  integer  variables  and  fixed  charges.  Unlike  the  binary  decision  variables 
in  GAP,  which  restricts  the  assignment  of  an  operation  to  only  one  machine  group,  partial 
assignments  of  operations  to  several  machine  groups  is  permitted  in  our  problem.  The 
machine  capacity  constraints  include  integer  decision  variables  to  determine  additional 


machine  requirements  in  order  to  meet  the  production  schedule.  Since  the  number  of 
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additional  machines  is  not  restricted,  a feasible  solution  to  the  problem  always  exists.  The 
fixed  costs  imposed  on  the  procurement  decision  variables  relate  the  problem  to  a Fixed 
Charge  Transportation  Problem  as  well.  Given  the  number  of  machines  in  each  machine 
group  and  capacity  consumption  of  each  operation  type  in  each  machine  group,  our 
objective  is  to  find  the  assignment  scheme  of  individual  operations  to  machine  groups 
that  minimizes  the  total  monthly  operating  cost  of  machines  processing  these  operations 
and  the  procurement  cost  of  additional  machines  capitalized  monthly.  To  the  best  of  our 
knowledge  no  research  exists  in  the  literature  that  addresses  this  type  of  problem  with 
integer  and  continuous  variables  involving  both  assignment  and  capacity  decisions.  Only, 
Askin,  Selim,  and  Vakharia  (1997)  approach  a similar  problem  in  the  operations 
assignments  phase  of  their  methodology  for  designing  flexible  cellular  systems  where  the 
resource  capacities  may  vary.  They  present  two  greedy  heuristics  to  find  the  resource 
capacity.  Since  they  do  not  permit  partial  assignments,  their  problem  reduces  then  to  the 
GAP. 

As  we  mentioned  earlier,  certain  operation  types  are  repeated  several  times  during 
the  manufacturing  process,  allowing  wafers  to  visit  the  same  machine  group  more  than 
once.  Since  the  specifications  of  the  process  performed  on  the  wafer  differ  from  one  visit 
to  another,  each  visit  is  referred  to  as  a distinct  step  and  each  processing  step  is  given  a 
specific  name.  This  name  indicates  both  the  wafer  type/s  and  the  process  to  be  performed 
on  the  wafer  type/s.  Throughout  this  chapter  wafer  type/s  at  a particular  processing  step 
is/are  referred  to  as  an  “operation”.  In  this  environment,  we  focus  on  the  assignment 
problem  of  operations  to  alternative  machine  groups  and  we  consider  a population  of 
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machines  that  are  only  capable  of  processing  these  operations.  Machines  have  different 
manufacturing  characteristics  according  to  the  process  and  product  involved.  While 
newer  machines  are  normally  capable  of  processing  advanced  products  as  well  as  older 
products,  older  machines  can  only  process  older  products,  may  require  longer  processing 
times  and  may  have  lower  yields  and  lower  utilization  levels  due  to  more  maintenance 
requirements.  Machines  having  the  same  characteristics  are  gathered  into  machine  groups 
and  each  machine  group  represents  a set  of  parallel  identical  machines  with  the  same 
processing  capabilities.  A machine  group  may  be  constituted  of  one  or  more  machines  in 
an  existing  manufacturing  system.  The  number  of  machines  is  zero  in  the  case  where  a 
new  machine  group  is  formed.  In  what  follows  a discussion  of  the  data  requirements  and 
mathematical  formulation  of  the  assignment  problem  is  presented. 

2.2.1.  Data  Requirements 

The  set  of  machine  groups,  1=  {1, . . . , m),  set  of  operations,  J-  {\, ...  ,n),  and 
the  number  of  machines  in  each  machine  group,  Nh  are  given.  For  an  existing 
manufacturing  system,  N,  is  the  number  of  machines  currently  available  in  each  machine 
group  i.  For  a new  machine  group  i,  Nt  is  zero. 

A matrix  [P]  is  defined  as  the  capability  matrix:  the  binary  entry  ptJ  is  1 if  machine 
group  i can  perform  operation  j and  0 otherwise.  The  operations  may  be  split  among 
machine  groups  and  can  be  performed  in  any  order.  Daily  scheduled  volume  of 
operations  (VJ),  the  yield  of  machine  type  i with  respect  to  operation  j ( Y j),  and  the 


30 


Table  2.1.  Illustrative  data  for  four  operations  processed  on  ion  implanters 


Machine  Type 

Quantity 

Available 

Net  Production 
Hrs/24  Hr 

Utilization 

Ion  Implanter  1 63 

3 

19.2 

0.80 

Operation 

Number 

Scheduled 

Volume 

Production 

Capacity 

% Yield 

Number  of 
machines  needed 

5 

95 

544 

98.12 

0.18 

11 

244 

1022 

95.26 

0.25 

16 

169 

898 

92.96 

0.20 

23 

44 

643 

96.37 

0.07 

quantity  of  operation  j (in  other  words,  number  of  wafers)  that  machine  group  i can 
process  per  day  (PCy)  are  assumed  to  be  known.  PC,j  is  computed  as  a function  of  the 
raw  processing  time,  setup  time,  and  operator  intervention  time.  Utilization  ([/,)  is  the  net 
production  hours  of  the  machine  group  i during  the  24  hour  period.  Given  these  data,  the 
number  of  machines  of  type  i required  to  process  operation  j ( cty)  is  computed  with  the 
following  formula,  assuming  that  all  machines  are  available  100%  of  the  time: 

v, 

“ I -P“ 


The  data  for  four  operations  are  displayed  in  Table  2.1  as  an  example.  Also  shown  is  a 
machine  group  consisting  of  three  ion  implanters  that  can  process  these  operations.  The 
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monthly  operating  cost  of  each  machine  of  type  i to  process  operation  j (cy)  and  monthly 
capitalized  procurement  cost  of  an  additional  machine  of  type  i ( F ’■)  are  also  known. 

2,2.2.  Problem  Formulation 

A comprehensive  mathematical  model  of  the  operation-machine  type  assignment 
problem  with  integer  and  continuous  decision  variables  and  fixed  charges  on  resources  is 
as  follows: 


(MIAP)  min  0) 

iel  jeJ  iel 

s.t.  Y,xij  =1’  JeJ  (2) 

iel 

X aiixii  £(Nt+  Qi )Ut,  iel  (3a) 

jeJ 

0 — xy  — Pij  > iel,  jeJ  (4) 

Q(  > 0 and  integer,  iel  (5) 


The  decision  variable  x y (0  <xy  < 1)  indicates  the  fraction  of  operation  j assigned  to 
machine  group  i and  integer  variable  Q{  represents  the  number  of  additional  machines 
required  in  machine  group  i.  atj  is  the  fraction  of  machine  type  i required  to  perform 
operation  j,  assuming  that  all  of  operation  j is  processed  on  machine  group  i.  Ut  can  be 
used  as  a parameter  to  control  the  amount  of  slack  capacity  built  into  the  system  (0  < Uj 
< 1).  For  a given  £/,  we  define  ay  = ay  / U,  and  replace  constraint  (3a)  with 
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IXxy  ~Q‘  ~N‘  ’ /6/  (3b) 

In  the  formulation  of  the  Mixed  Integer  Assignment  Problem  (MIAP),  the 
objective  function  (1)  minimizes  total  monthly  operating  cost  of  machine  groups  and 
monthly  procurement  cost  of  additional  machines.  The  first  constraint  set  (2)  enforces  all 
scheduled  volume  for  an  operation  to  be  completely  assigned.  Constraint  set  (3b)  ensures 
that  there  is  sufficient  machine  capacity  in  each  machine  group  to  process  the  assigned 
operations  considering  machine  utilization  and  machine  duplication.  Constraints  (4) 
prevent  assignments  of  operations  to  machine  groups  that  are  not  capable  of  processing 
them  and  constraints  (5)  impose  integrality  restriction  on  the  number  of  additional 
machines.  Hence,  this  model  finds  the  optimum  assignment  scheme  of  operations  to 
machine  groups  to  minimize  total  operating  and  procurement  costs  of  machines. 

In  the  next  section,  we  apply  Lagrangian  relaxation  and  Lagrangian 
decomposition  techniques  with  subgradient  optimization  to  MIAP  to  obtain  lower  bounds 
on  the  optimal  solution.  Lagrangian  relaxation  offers  an  alternative  to  linear  programming 
relaxation  for  finding  good  bounds  on  the  optimal  solutions.  Lagrangian  decomposition  is 
an  extension  to  Lagrangian  relaxation  approach  and  may  yield  improved  bounds  on  the 
optimal  solutions.  An  efficient  solution  procedure  is  also  developed  to  find  a good 


assignment  scheme  for  the  problem. 
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2.3.  Lower  Bounding  Methods 

Lagrangian  relaxation  and  Lagrangian  decomposition  techniques  are  discussed  in 
this  section  to  obtain  good  lower  bounds  on  the  optimal  solution  to  MIAP.  These  two 
relaxation  approaches  have  been  successfully  employed  by  numerous  researchers  in  the 
past  in  many  integer/mixed  integer  programming  applications. 

2.3.1.  Lagrangian  Relaxation 

Lagrangian  relaxation  is  a well  known  way  to  circumvent  the  direct  consideration 
of  “complicating”  constraints  in  an  optimization  model  and  has  been  widely  used  in  large 
scale  mathematical  programming.  Lagrangian  relaxation  is  formed  by  multiplying  the 
complicating  constraints  with  corresponding  penalty  costs  (Lagrangian  multipliers)  and 
including  them  into  the  objective  function.  By  dualizing  those  constraints,  we  obtain  a 
problem  that  is  easy  to  solve  and  whose  optimal  value  gives  a lower  bound  on  the  optimal 
value  of  the  minimization  problem.  Although  the  method  does  not  guarantee  that  a 
solution  will  be  found  for  all  problems,  it  is  often  worth  trying  since  it  is  so  simple 
compared  to  other  methods.  The  reader  is  referred  to  Geoffrion  (1974)  and  Fisher  (1981) 
for  a detailed  and  excellent  study  on  the  Lagrangian  relaxation  technique.  Cattrysse  and 
Van  Wassenhove  (1992)  provide  an  extensive  survey  of  several  Lagrangian  relaxation 
methods  developed  in  the  literature  to  solve  the  GAP. 
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The  relaxation  to  our  assignment  problem  is  obtained  by  dualizing  constraint  set  (2).  The 
yielding  Lagrangian  relaxation  LR(/l)  with  Lagrangian  multiplier  vector  X has  the 
objective  function 


where  X is  unrestricted  in  sign.  This  problem  nicely  decomposes  into  m continuous 
knapsack  subproblems,  one  for  each  machine  group  i.  A procedure  to  solve  this  set  of 
problems  is  as  follows: 

Step  1:  Order  xp’s  in  the  increasing  order  of  ( Cy  ay  - Xj)  / ay  value. 

Step  2:  Proceed  along  the  list:  if  ( Cy  a y - Xj  ) I cij  is  positive  or  zero,  the 


corresponding  Xy  is  set  equal  to  0.  If  ( Cy  ay  - Xj)  I ay  is  negative, 
the  corresponding  Xy  is  set  to  equal  1 assigning  the  whole 
production  volume  to  machine  group  i,  provided  there  is  enough 
capacity  left  at  that  machine  group,  taking  into  account  the  other 
operations  who  have  already  been  assigned  to  machine  group  i. 
However,  if  there  is  not  sufficient  capacity  at  machine  group  i to 
assign  operation  j completely,  set  xy  as  high  as  possible  to  exhaust 
capacity  at  machine  group  i. 


Step  3:  Compare  the  trade-off  between  the  cost  of  adding  an  additional 


machine  capacity  and  assigning  remaining  operations  following 
the  procedure  in  Step  2.  If  the  difference  is  negative,  increase  Qi  by 
1 and  assign  operations  as  described  in  Step  2. 


minimize 


+ I ^0,  + E * 


iel 
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Step  4:  Repeat  Step  3 incrementing  Qj  until  the  objective  function  value 

ceases  to  decrease. 

The  preliminary  analysis  of  this  relaxation  showed  that  the  structure  of  the 
objective  function  coefficients  largely  affects  the  optimal  value  and  weak  lower  bounds 
may  be  attained  in  some  instances.  Therefore,  we  considered  adding  a lower  bound  on  the 
total  machine  requirements  over  all  machine  groups  and  dualizing  it  along  with  the 
capacity  constraints.  The  new  constraint  is  as  follows: 

(6) 

iel 

We  observed  that  this  relaxation  method  could  substantially  improve  the  lower  bound. 
The  integer  scalar  K can  be  obtained  using  the  following  two  methods: 

• Solve  the  linear  programming  problem  that  minimizes  total  number  of 
additional  machines  required  subject  to  (2)-(4)  and  round  up  the 
objective  function  value. 

• Consider  the  best  case  scenario:  assign  each  operation  j to  the  machine 
group  i where  it  consumes  least  amount  of  resources  regardless  of  the 
machine  group’s  capacity  but  subject  to  its  capability.  Identify  the 
number  of  additional  machines  required  in  each  machine  group  i.  Round 
up  the  difference  between  the  total  number  of  machines  needed  and  total 
number  of  machines  currently  available. 
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Lagrangian  relaxation  LR(/l,/v),  X unrestricted  and  juz  0,  with  respect  to  constraints  (2) 
and  (6)  is  as  follows: 


LRyW/y)  min 

HH(Cua0  ~Xi)xij 

iel  jeJ 

+ Z + ’£*,  + km 

iel  jeJ 

(7) 

S.t. 

iel 

Xy-Q^Ni, 

i G / 

(8) 

0<Xy  <Py  , 

ie  I,  j ej 

(9) 

Qi  > 0 and  integer. 

i e I 

(10) 

We  use  a subgradient  optimization  procedure  to  update  the  Lagrangian 
multipliers.  Given  initial  values  of  X°  and  ju°,  we  generate  a sequence  of  Lagrangian 
multipliers  X1  and  /u 1 through  the  addition  of  a direction  vector  which  is  multiplied  by  a 
step  size  9'  where  9 ' is  a positive  scalar.  Hence,  the  updating  of  Lagrangian  multipliers  of 
the  constraint  set  (2)  is  done  according  to  the  equation 

f \ 


t-\ 


+ 9‘ 


'-I 


j G J 


V iel  y 


and  Lagrangian  multiplier  of  the  constraint  (6)  is  updated  according  to  the  equation 


H ' = max  i 0 ;//'*'+  9 1 


k- I a 

V is/ 


The  step  size  9‘  is  updated  at  each  iteration  t using  the  following  equation: 
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6'  = S' 


(UB  -LB  (AM ; 


r 

f \ 

2 / 

ji 

1 - Z xu 

+ 

/ JeJ 

\ ‘el  7 

V 

K - z 0, 


ie/  / 


The  step  size  0‘  depends  on  the  parameter  5‘  { 0 < 5‘  <2  ),  on  the  gap  between  the 
current  lower  bound  LB(TM,  /u M)  and  the  estimated  minimum  value  of  the  objective 
function  of  the  relaxed  problem,  which  is  approximated  by  the  upper  bound  UB  obtained 
by  applying  a heuristic  method,  and  on  the  Euclidean  norm  of  the  deviations  in  the 

relaxed  constraints  (2)  and  (6).  The  sequence  S'  is  determined  by  setting  c>°  = 2 and  by 
halving  8l  whenever  /u  M)  does  not  increase  after  a fixed  number  of  iterations. 

We  terminate  this  procedure  when  one  of  the  following  stopping  criteria  is  met: 

(1)  An  iteration  number  limit, 

(2)  Maximum  gap  between  the  lower  and  upper  bounds, 

(3)  A limit  on  the  value  of  the  Euclidean  norm  of  the  deviations. 

2.3.2.  Lagrangian  Decomposition 

The  Lagrangian  decomposition  approach  was  introduced  by  Jomsten  and  Nasberg 
(1986).  The  idea  behind  is  that  constraints  sets  (2)  and  (3b)  are  totally  different,  but  both 
sets  of  constraints  have  a special  structure.  The  traditional  Lagrangian  relaxation  dualizes 
one  of  the  constraint  set  and  keeps  the  other.  However,  we  may  be  able  to  obtain  better 
lower  bounds  on  the  optimal  objective  function  value  by  keeping  both  constraint  sets  in 
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the  formulation.  Consider  the  following  reformulation  of  MIAP  with  the  lower  bounding 
restriction  K on  the  sum  of  the  total  number  of  additional  machines  required: 


(AMIAP)  min 
s.t. 


P YLcijaijxii  +r 

iel  jeJ 

iel  jeJ  iel 

(11) 

ii 

j eJ 

(12) 

I ]a0x0  -Qi^Ni  » 

iel 

i e / 

(13) 

Za**. 

16/ 

(14) 

xa=y u> 

iel,  jeJ 

(15) 

xg.y9  ^ Py  > 

iel,  j eJ 

(16) 

Qi  > 0 and  integer, 

i e I 

(17) 

where  (3  and  / are  nonnegative  parameters.  MIAP  and  AMIAP  are  equivalent  since  the 
objective  function  values  are  related  in  the  following  way: 

v(MIAP)  = (/?+/)  v( AMIAP) 

where  v(P)  is  the  objective  function  value  of  problem  (P).  Notice  that  v(MIAP)  = 
v( AMIAP)  when  J3+  y=\.  When  the  constraint  sets  (14)  and  (15)  are  relaxed, 
Lagrangian  relaxation  separates  nicely  into  two  problems,  one  problem  with  the  x- 
variables  only  and  another  with  /-variables  only: 
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]JDx(A,/u) 

min 

YL(P  cijaij-*-ij)xij 

iel  jeJ 

iel 

(18) 

s.t. 

iel 

i e I 

(19) 

0<Xy  < Py  , 

iel,  j eJ 

(20) 

Qi  > 0 and  integer, 

i e I 

(21) 

LD y(A) 

min 

iel  jeJ 

(22) 

s.t. 

II 

W5 

jeJ 

(23) 

0iy9sp,, 

iel,  j eJ 

(24) 

where  //  is  the  nonnegative  Lagrangian  multiplier  vector  for  constraint  (14)  and  A, 
unrestricted,  is  the  Lagrangian  multiplier  vector  for  constraint  set  (15).  Then,  the 
objective  function  value  of  AMIAP  is 

v(AMIAP)  = v{LDx{A,/li))  + v(LD y(A))  + K/u 
\JDx(A,/j)  decomposes  into  m continuous  knapsack  problems,  LDjc ,{A,ju),  which  can  be 
solved  using  the  procedure  described  in  the  previous  section.  LDy(/l),  on  the  other  hand, 
can  be  separated  into  n trivial  generalized  upper  bound  problems,  LDv/T),  where  the 
assignment  decision  variables  are  naturally  binary.  The  solution  procedure  is  as  follows: 
For  all  j e J,  determine  argmin  i€j  {ycy  atj  + Ay  } and  set  the  associated  ytj  = 1 . 
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Remaining  y^s  are  set  to  zero.  Lagrangian  multipliers  in  both  problems  are  updated  using 
the  subgradient  optimization  procedure  as  described  in  Lagrangian  relaxation  method. 


2.4,  A Solution  Procedure  for  Operation-Machine  Type  Assignment  Problem 


In  this  section,  we  propose  a four  phase  solution  procedure  to  the  Operation- 
Machine  Type  Assignment  Problem.  An  outline  of  the  procedure  is  given  in  Figure  1.  In 
Phase  I we  investigate  linear  programming  relaxation  of  the  problem.  If  the  variables 
corresponding  to  machine  duplication  decisions  are  zero  or  integer  then  LP  solution  is 
optimal.  Since  this  is  very  rarely  the  case  we  use  the  LP  relaxation  to  obtain  a lower 
bound  on  the  objective  function  and  an  initial  solution  for  the  improvement  procedure. 
Furthermore,  if  the  problem  is  feasible  with  the  number  of  additional  machines  set  equal 
to  zero,  the  solution  gives  an  upper  bound  on  the  value  of  the  objective  function. 

7 


PHASE  I: 

PHASE  II: 
PHASE  III: 

PHASE  IV: 


Solve  LP  relaxation  of  MIAP  (LP1).  If  the  solution  to  LP1  is 
feasible  to  MIAP  then  it  is  also  the  optimal  solution,  stop. 
Otherwise,  LP1  solution  gives  a lower  bound  to  v(MIAP). 

If  linear  programming  problem  (LP2)  where  (9,  = 0 for  Vi  el 
is  feasible  then  the  optimal  solution  value  is  an  upper  bound 
to  v(MIAP). 

Use  a greedy  heuristic  to  reduce  the  cost 
Use  Lagrangian  relaxation  (LR)  and/or  Lagrangian 
decomposition  (LD)  method  to  improve  the  lower  bound. 
Solve  linear  programming  problem  for  fixed  values  of  Q{ 
obtained  in  Phase  II  to  find  the  optimal  assignment  scheme 
and  update  the  upper  bound. 


Figure  2. 1 . Operation-Machine  Type  Assignment  Methodology 
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Phase  II  is  a greedy  heuristic  consisting  of  an  improvement  procedure  where  we 
consider  pair  wise  exchanges.  Starting  from  an  initial  feasible  solution,  the  procedure 
evaluates  all  possible  exchanges  between  pairs  of  machines  and  makes  an  exchange  if  it 
improves  the  value  of  the  objective  function.  We  select  the  exchange  providing  the 
maximum  cost  improvement.  The  initial  solution  is  obtained  by  solving  the  LP  relaxation 
of  the  problem  and  rounding  up  the  decision  variables  corresponding  to  additional 
machine  requirements  to  meet  the  integrality  condition  of  the  original  problem.  The 
objective  function  value  updated  according  to  integer  values  gives  an  upper  bound.  Since 
the  number  of  machines  are  rounded  up,  there  exists  at  least  one  machine  group  with 
slack  capacity.  Therefore,  in  the  next  step  we  investigate  the  pairs  of  machine  groups  to 
check  whether  removing  an  operation,  partially  or  entirely,  from  one  machine  group  and 
reassigning  it  to  another  to  reduce  the  cost.  This  may  yield  a reduction 

• in  the  operating  cost,  if  it  is  cheaper  to  process  that  operation  in  the  second 
machine  group, 

• in  the  procurement  cost,  if  the  removal  of  that  operation  results  in  a saving  of  one 
machine, 

• both  in  operating  and  procurement  costs. 

After  all  pairs  of  machine  groups  are  investigated,  the  pair  that  achieves  the 
maximum  cost  reduction  is  selected  and  the  assignments  and  machine  requirements  are 
updated  accordingly.  Note  here  that  capacity  consumption  of  a certain  operation  may 
differ  from  one  machine  group  to  another  since  machines  have  different  processing 
capabilities.  Therefore,  capacity  consumption  comparisons  between  machine  groups  are 
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performed  accordingly,  so  are  the  machine  requirements  in  each  group  updated.  Table  2.2 
depicts  an  example  where  0.66  and  0.12  units  of  additional  machines  are  needed  in 
machine  group  i and  machine  group  j,  respectively  (i.e.  1 machine  is  needed  in  each 
group).  If  a fraction  of  operation  j is  transferred  at  no  additional  procurement  cost  from 
machine  group  / to  machine  group  i where  the  operating  cost  is  lower  then  an 
improvement  on  the  total  cost  is  possible.  In  addition,  if  that  fraction  of  operation  j 
requires  0.12  units  of  machine  capacity  /,  total  cost  is  further  improved  by  saving  one 
machine.  In  that  case,  the  additional  machine  requirement  reduces  to  0 and  total 
procurement  cost  decreases  by  the  cost  of  one  machine  of  type  /.  Table  2.3  illustrates  this 
case  where  25%  of  operation  j is  removed  from  machine  group  / and  reassigned  to 
machine  group  i. 


Table  2.2.  Capacity  allocation  prior  to  operation  exchange 


Operation  j . 

Operation  k 

Total 

Available 

machines 

Machines 

required 

Mach.  Group  i 

0.18 

0.86 

. 2.66 

2 

1 

Mach.  Group  / 

0.82 

0.00 

. 3.12 

3 

1 

Table  2.3.  Capacity  allocation  after  the  operation  exchange 


Operation  j . 

Operation  k 

Total 

Available 

machines 

Machines 

required 

Mach.  Group  i 

0.43 

0.86 

. 2.82 

2 

1 

Mach.  Group  / 

0.57 

0.00 

. 3.00 

3 

0 
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At  the  end  of  Phase  II,  the  final  values  of  the  additional  machine  requirements  are 
obtained  to  build  the  new  machine  plan.  In  Phase  III,  we  try  to  improve  the  lower  bound 
using  Lagrangian  relaxation  and/or  Lagrangian  decomposition  method  as  described  in 
Section  3.  Finally,  in  Phase  IV  we  solve  a linear  programming  problem  to  find  the 
optimal  assignment  scheme  for  the  achieved  machine  plan  and  we  update  the  upper 
bound  on  the  optimal  value  of  the  objective  function. 

In  what  follows  is  a detailed  description  of  cost  reduction  algorithm  in  Phase  III. 
We  first  introduce  the  following  notation  for  the  data  and  variables  employed  in  the 
algorithm: 

qt  number  of  additional  type  i machines  required  (may  be  decimal) 

Diff  reduction  in  Xy 

MTj  number  of  type  i machines  to  be  transferred 

MNj  number  of  total  type  i machines  needed  (may  be  decimal) 

MEk  excess  capacity  of  type  k machines 


MA  k - MN  k 
floor  (MN  k)-MN  k. 


otherwise 


if  MAk  > MN  k 


MAi 


number  of  available  type  i machines 


TC 


total  operating  and  procurement  costs 


ATC 


change  in  total  cost 


AOCjkj  change  in  operational  cost  by  transferring  operation  j,  partially  or 


entirely,  from  machine  type  i to  machine  type  k 
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ATCikj 


DiffBUY 

ATCbuy 

DiffNOTBUY 

ATCNOtbuy 

ceil(x) 

floor(x) 


change  in  total  cost  by  transferring  operation  j,  partially  or  entirely, 
from  machine  type  i to  machine  type  k 
reduction  on  xy  if  a new  machine  of  type  k is  added 
change  in  total  cost  if  a new  machine  of  type  k is  added 
reduction  on  xy  if  no  machine  is  added  to  machine  group  k 
change  in  total  cost  if  no  machine  is  added  to  machine  group  k 
smallest  integer  number  greater  than  or  equal  to  x 
greatest  integer  number  less  than  or  equal  to  x 


Greedy  Heuristic: 

Step  0:  Round  up  qt  obtained  in  LP1  to  integer  values  to  obtain  Qt : Q,  = ceil(#,). 

Update  the  objective  function  value  TC. 

Step  1:  Let  / = { i e I\  Q{ > 0 }.  If/  = 0,  go  to  Phase  IV. 

For  Vi  el , go  to  Step  2. 

Step  2:  Set  MTt  = MNt  - floor(MV,).  If  MTt  = 0,  let  MTt  = 1 . 

Go  to  Step  3. 

Step  3:  Let  J = {j  e J\ Xy  > 0 } 

For  \/k  e I-  {/}  and  V/  e J repeat  Step  4-6. 

Step  4:  Set  Axy  = MTt  / ay  and  A VCikj  = ckj  - Cy. 


Go  to  Step  5. 
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Step  5:  If  ( Axjj  < xtj ) set  Drop  = 1 . Otherwise,  set  Drop  = 0 and  Axy  = x^. 

Compute  MEk. 

Step  6:  If  ( Ax^  < MEk  / akj),  let  Diffikj  = AxtJ,  and 

ATCikj  = (AVCikj  * Diffikj)  - {Fi  * Drop) 

Otherwise,  let 

Diff buy  = Axy  < ATC buy  = ( AVCikj  * Diff Buy)  - (Fi  * Drop)  + Fk 

Diff notbuy  = MEk  / akj , ATC notbuy  - (AV Cikj  * Diff  notbuy) 

ATCtkj  = min  { ATC  buy  > ATC  notbuy  } 

\ Delta  BUY  , if  ATC  — ATC  BUY 
Let  Diff  =< 

1 7 [ Delta  NOTBUY  , otherwise 

Step  7:  Find  /*,/,  k*  = argmin  ATCikj . 

Let  ATC  = ATCi*k*j*  and  Diff  = Diff*k*j * . 

Step  8:  If  ATC  = 0,  stop,  no  improvement  on  TC  is  possible.  Otherwise: 

TC  =TC+  ATC* 

Xj*j*  x i *j  * - Diff,  xk*j*  xk*j*  * Diff 

MN, * = MNj*  - Diff . aj*j*,  MNk*  = MNk*  - Diff.  ak*j* 

If  (Diff  = MTj*/  at*j*)  qt*  = qt*  - 1 

If  (Diff > MEkJ  akV* ) qk*  = qk*  + 1 

Let  Qi  = ceil(<7,)  and  Qk  = ceil(^).  Go  to  Step  1 . 
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2.5.  Experimental  Analysis 


The  Operation-Machine  Type  Assignment  Methodology  is  coded  in  C 
programming  language.  CPLEX®  optimization  software  was  used  to  solve  LP  problems. 
A series  of  computational  results  was  carried  out  on  a PC  with  300  MHz  Pentium  II 
processor  using  several  sets  of  randomly  generated  problems  and/or  problems  based  on 
industry.  Since  the  industry  data  do  not  include  operating  and  procurement  cost  figures, 
those  costs  are  generated  randomly.  Considering  different  sized  problems,  the 
performance  is  evaluated  using  various  combinations  of  cost  structures.  The  procurement 
cost  coefficient  Fj  is  generated  from  uniform  distributions:  D,(F/)  ~ U(10,  12), 
D2(F/)~U(10,  15),D3(F,-)~U(10,20),  and  D4(F;)~U(10,  40).  The  operating  cost 
cij  is  a percentage  of  the  procurement  cost  and  is  drawn  from  uniform  distributions: 
D,(c,y)  ~ U(0.1,  0.3),  D2(cij)  ~ U(0.1,  0.9),  and  D 3(c//)  ~ U(0.1,  1.9).  In  other  words,  the 


Table  2.4.  Summary  of  results 


Probl 

Prob2 

Prob3 

Prob4 

Prob5 

Prob6 

Prob7a 

(Z*-LB)/LB 

Best 

Average 

Worst 

0.00% 

0.20% 

4.45% 

0.00% 

0.74% 

5.72% 

0.00% 

1.09% 

10.60% 

0.00% 

0.50% 

3.07% 

0.00% 

0.39% 

2.85% 

0.00% 

0.22% 

2.13% 

0.01% 

1.28% 

7.25% 

Best 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

(UB-Z*)/Z* 

Average 

0.54% 

1.90% 

1.74% 

1.48% 

1.34% 

1 .02% 

Worst 

7.44% 

8.81% 

9.54% 

8.65% 

7.13% 

5.52% 

a Z*  could  not  been  found  for  Prob7,  the  results  are  based  on  (UB-LB)/LB 
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operating  cost  is  on  the  average  20%  of  the  procurement  cost  using  D,(cy),  50%  using 
D,(cy),  and  100%  using  D,(cjj).  We  consider  three  capacity  utilization  ratios:  0.60,  0.75, 
and  0.90. 

We  tested  the  solution  procedure  on  seven  types  of  problems  over  a wide  range  of 
values  for  the  number  of  machine  groups  and  number  of  operations.  Each  problem  type 
was  solved  five  times  with  the  same  input  structure  to  achieve  a reasonable  level  of 
confidence  to  validate  the  results  (4  x 3 x 3 = 36  input  parameters,  36  x 5 = 180  runs  for 
each  problem  type).  A summary  of  the  results  is  presented  in  Table  2.4.  Optimal  solution 
for  problem  type  7 cannot  be  obtained  because  of  the  size  of  the  problem  and  computer 
memory  limitation.  For  all  other  problems,  the  optimal  solution  was  achieved  at  least 
once  using  Lagrangian  relaxation/decomposition  technique  and  the  greedy  heuristic.  The 
gap  between  the  optimal  solution  and  lower  bound  is  less  than  1.1%  on  the  average  in  all 
problem  instances  and  10.6%  in  the  worst  case;  and  the  gap  between  the  upper  bound  and 
optimal  solution  is  less  than  1.9%  on  the  average  in  all  problem  instances  and  9.54%  in 
the  worst  case. 

The  results  show  that  better  lower  bounds  and  upper  bounds  are  achieved  when 
the  number  of  operations  in  the  problem  is  greater  than  the  number  of  machine  groups. 
Figure  2.2  and  2.3  illustrate  the  variation  in  mean  gaps  with  respect  to  procurement  costs 
and  operating  costs,  respectively,  at  different  utilization  ratios.  A pattern  cannot  be 
apparently  observed  between  those  parameters.  On  the  other  hand,  although  the  gap 
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Means  for  Utilization  = 0.90 


4.000% 

N 3.000% 

S'  2.000% 

CQ 

5 1.000% 

0.000% 

( 10,12)  ( 10,15)  (10,20)  (10,40) 
Procurement  Cost 

( 0.1, 0.3)  ( 0.1,  0.9 ) ^-  ( 0. 1,  1 .9) 


Figure  2.2.  Mean  gaps  over  procurement  costs  for  different  levels  of  utilization 
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Means  for  Utilization  = 0.60 


o 0.800% 

^ 0.600% 

3 0.400% 

*•  0.200% 

& 0.000% 

(0.1, 0.3)  (0.1,  0.9)  (0.1,  1.9) 

Operating  Cost 


Means  for  Utilization  = 0.60 

» 2.000% 

1.500% 

N 1.000% 

m 0.500% 

- 0.000%-- i 

( 0.1, 0.3  ) (0.1, 0.9)  (0.1,  1.9) 


Operating  Cost 


-*-(  10,  12)  -m-(  10,  15) 
-*-(  10,20)  -*-(  10,40) 


-♦-(  10,  12)  -®-(  10,  15) 
-*-(  10,20)  -*-(  10,40) 


Figure  2.3.  Mean  gaps  over  operating  costs  for  different  levels  of  utilization 
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Means  for  Operating  Cost  ~ U(0. 1,1.9) 


-*-(10,12)  ( 10,  15  ) 

— ( 10,  20  ) -*-(  10,40) 


Means  for  Operating  Cost  ~ U(0. 1,1.9) 

2.000%  , 

N 1.500% 

S'  1.000% 

i 

§ 0.500% 

0.000%  | 

0.60  0.75  0.90 

Procurement  Cost 


(10,12) 
• ( 10,  20) 


-(  10,  15) 
-(  10,  40) 


Figure  2.4.  Mean  gaps  over  utilizations  for  different  procurement  and  operating  costs 
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between  the  optimal  solution  and  lower  bound  seems  to  increase  as  utilization  increases 
we  cannot  make  any  statement  for  the  gap  between  the  upper  bound  and  optimal  solution. 
Figure  2.4  shows  the  change  in  mean  gaps  with  respect  to  utilization  ratios.  Detailed 
results  of  all  seven  problems  are  given  in  Appendix  A. 


CHAPTER  3 

MULTI-PERIOD  CAPACITY  ALLOCATION  WITH  MACHINE  DUPLICATION  IN 
SEMICONDUCTOR  MANUFACTURING 

3.1.  Introduction 

Most  of  the  recent  research  in  semiconductor  manufacturing  has  focused  on 
demand  satisfaction,  maximization  of  equipment  utilization,  minimization  of  production 
costs,  and  maximization  of  throughput  with  some  capacity  constraints.  Linear 
Programming  is  often  proposed  as  a tool  for  production  planning  and  scheduling  in 
semiconductor  manufacturing.  However,  LP  formulations  can  be  very  large  for  large 
organizations  with  complex  production  environments  such  as  semiconductor  industry 
requiring  a very  long  time  to  generate  the  input  data  files  that  need  to  be  fed  into 
mathematical  planning  software  and  huge  amount  of  memory  and  disk  space  to  store 
these  data.  Golovin  (1986)  points  out  the  difficulty  of  selecting  an  appropriate  objective 
function  in  the  planning  of  the  semiconductor  manufacturing  and  shows  that  a detailed 
formulation  of  the  problem  in  an  integrated  manner  based  upon  a mathematical 
programming  model  is  intractable  and  requires  data  which  cannot  be  obtained  reliably. 
He  proposes  a hierarchical  approach  such  as  that  suggested  by  Bitran  et  al.  (1981,  1982). 

Leachman  (1993)  gives  an  optimization-based  corporate-level  production 
planning  system  that  includes  multiple  facilities  and  treats  entire  production  processes  in 
each  facility  as  integral  entities.  His  model  generates  capacity-feasible  start  and  out 
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schedules  for  each  manufacturing  facility  in  the  company.  Leachman  and  Carmon  (1992) 
analyze  capacity  of  semiconductor  production  facilities  in  which  manufacturing 
operations  may  be  performed  by  alternative  machine  types.  This  type  of  operations  can  be 
performed  by  different  machine  types  with  different  processing  times.  They  propose  a 
modeling  technique  that  greatly  reduces  the  size  of  the  LP  problems.  Hung  and  Wang 
(1997)  present  an  alternative  model  for  formulating  the  alternative  machine  capacitated 
planning  problem  as  an  extension  of  the  technique  proposed  by  Leachman  and  Carmon 
(1992).  Their  technique  considerably  reduces  the  size  of  the  LP  problem  for  bin 
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Figure  3.1.  Simple  silicon  TTL  Integrated  Circuit  process  flowchart 
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allocation  planning  in  semiconductor  manufacturing,  thus  saving  substantial  amount  of 
computer  resources  required.  The  interested  reader  is  referred  to  Uzsoy  et  al.  (1992)  for  a 
comprehensive  review  of  the  research  in  the  semiconductor  production  planning  and 
scheduling. 

The  alternative  machine  type  capacity  allocation  problem  with  duplicated 
machines  in  a complex  semiconductor  manufacturing  environment  is  a difficult  task. 
Machines  have  different  manufacturing  characteristics  according  to  the  process  and 
product  involved.  While  newer  machines  are  normally  capable  of  processing  advanced 
products  as  well  as  older  products,  older  machines  can  only  process  older  products,  may 
require  longer  processing  times  and  may  have  lower  yields  and  lower  utilization  levels 
due  to  increasing  maintenance  requirements.  In  this  chapter,  we  address  this  problem  by 
recognizing  capacity  limitations  of  the  individual  machines  as  well  as  reducing  operating 
and  investment  costs  related  to  the  machines  and  inventory  holding  costs.  We  investigate 
the  tradeoff  between  allocating  a constant  capacity  for  the  planning  horizon  and  holding 
inventory  and  propose  a Lagrangian  relaxation  based  heuristic  approach  to  solve  the 
problem. 

The  remainder  of  this  chapter  is  organized  as  follows.  In  Section  3.2,  we  describe 
the  problem  and  formulate  a mathematical  model.  In  Section  3.3,  we  present  a 
Lagrangian-based  solution  procedure  to  obtain  lower  bound  and  a good  feasible  solution 
to  the  problem.  Section  3.4  contains  experimental  analysis  of  the  proposed  solution 
procedure  on  randomly  generated  data. 
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3.2.  Problem  Statement  and  Formulation 

Manufacturing  a computer  chip  is  a complex  process  involving  hundreds  of  steps 
and  requiring  from  a few  days  to  up  to  three  months  of  processing  time  to  complete.  The 
semiconductor  manufacturing  process  consists  of  various  operation  types  such  as 
oxidation,  deposition,  lithography,  etching,  ion  implantation,  etc.  performed  on  different 
batches  of  products  (wafers).  Some  of  these  operation  types  are  repeated  several  times  to 
build  different  layers  on  top  of  the  wafer  and  wafers  may  require  up  to  400  operation 
steps  at  a number  of  work  centers.  Each  wafer  type  follows  a particular  sequence 
throughout  its  fabrication  process  visiting  different  work  centers.  Specifically,  each  wafer 
makes  multiple  visits  to  the  same  work  center  at  different  points  in  the  fabrication 
process.  Such  a product  flow  is  known  as  a reentrant  product  flow  and  is  illustrated  in 
Figure  3.1.  The  common  practice  in  the  layout  of  work  centers  for  different  types  of 
processes  is  to  group  similar  operation  types  together:  furnaces  are  grouped  in  one  area, 
ion  implanters  in  another,  and  so  forth.  This  layout  requires  wafers  moving  back  and  forth 
between  work  centers  but  it  allows  the  utilization  of  the  same  equipment  to  process 
wafers  at  different  steps.  For  instance,  all  oxidation  processes  may  be  completed  on  the 
same  machine  or  in  the  same  work  center  although  the  wafers  are  required  to  travel  to 
other  work  centers  between  these  oxidation  process  steps. 

Since  the  specifications  of  the  process  performed  on  the  wafer  differ  from  one 
visit  to  another,  each  visit  is  referred  to  as  a distinct  step  and  each  processing  step  is  given 
a specific  name.  This  name  indicates  both  the  wafer  type/s  and  the  process  to  be 
performed  on  the  wafer  type/s.  Throughout  this  chapter  wafer  type/s  at  a particular 
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processing  step  is/are  referred  to  as  an  “operation.”  In  this  environment,  we  focus  on  the 
loading  of  operations  to  alternative  machine  types  where  machine  duplication  is  allowed 
for  certain  machine  groups.  We  consider  a population  of  machines  that  are  capable  of 
processing  these  operations  only.  Machines  having  same  characteristics  are  gathered  into 
machine  groups  and  each  machine  group  represents  a set  of  parallel  identical  machines 
with  the  same  processing  capabilities.  In  this  context,  machines  groups  are  classified  as 
primary  machine  groups  and  secondary  machine  groups  where  primary  machines  are  the 
most  efficient  machines  to  process  the  specified  operations  and  secondary  machines  are 
the  alternative  machines  in  case  additional  machine  capacity  is  needed.  Only  primary 
machines  may  be  procured  if  further  capacity  is  required  to  meet  the  production 
schedules. 

The  problem  of  modeling  the  capacity  of  alternative,  non-identical  resources 
arises  in  almost  all  manufacturing  environments  where  process  technology  is  evolving. 
Suprisingly,  published  research  on  this  topic  is  scant  (Leachman  and  Carmon,  1992).  In 
this  section,  we  develop  a mixed-integer  programming  formulation  to  minimize  operating 
cost  of  machines,  procurement  cost  of  primary  machines,  and  inventory  holding  cost.  For 
each  operation,  there  are  allocation  variables  to  spread  the  product-operation  production 
quantities  among  alternative  machine  groups.  Inventory  variables  and  integer  variables 
for  additional  machine  requirements  are  used  to  meet  the  production  volumes  at  each 
period.  Inventory  balance  constraints  are  formulated  to  guarantee  consistency  of 
production  volumes  between  consecutive  periods  and  resource  constraints  are  formulated 
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Table  3.1.  Illustrative  data  for  five  operations  processed  on  lithography  machines 


Machine  Type 

Quantity 

Available 

Net  Production 
Hrs/24  Hr 

Utilization 

PL  - 100 

5 

18.0 

0.75 

Operation 

Number 

Scheduled 

Volume 

Production 

Capacity 

% Yield 

Number  of 
machines  needed 

8 

212 

880 

97.44 

0.23 

17 

197 

1198 

96.38 

0.16 

24 

25 

476 

91.27 

0.05 

28 

52 

549 

93.44 

0.09 

39 

155 

1020 

95.71 

0.15 

both  to  satisfy  capacity  requirements  and  to  determine  additional  machine  requirements. 
Given  the  number  of  machines  in  each  machine  group  and  capacity  consumption  of  each 
operation  type  in  each  machine  group,  our  objective  is  to  find  the  optimal  assignment 
scheme  of  operations  to  machine  groups.  To  the  best  of  our  knowledge  no  research  exists 
in  the  literature  that  addresses  this  type  of  problem  with  integer  and  continuous  variables 
involving  both  assignment  and  capacity  decisions. 

Table  3.1  illustrates  sample  data  for  a machine  group  consisting  of  ion  implanters 
and  four  operations  that  can  be  processed  on  this  machine  group.  For  simplicity  of 
exposition,  we  will  assume  100%  yields  and  zero  initial  inventories.  We  introduce  the 
following  notation  for  the  data  concerning  our  capacity  allocation  problem: 
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t e T is  the  time  period  index,  where  t = 1, . . r, 
i e I is  the  machine  group  index,  where  i = 1, . . .,  m, 
j e J is  the  operation  index,  where  j = 1, . . .,  n, 

Nti  denotes  the  number  of  machines  currently  available  in  each  machine 
group  i at  period  t, 

[Z3]  is  the  capability  matrix:  the  binary  entry  is  1 if  machine  group  i can 
perform  operation  j and  0 otherwise, 

Vtj  is  the  daily  scheduled  production  volume  of  operation  j at  period  t, 

Ytij  is  the  yield  of  machine  type  i with  respect  to  operation  j at  period  t, 

PCtij  is  the  number  of  wafers  of  type  j that  machine  group  i can  process  daily 
at  period  t, 

U>i  denotes  the  utilization  of  the  machine  group  i at  period  t, 

( = Net  production  hours  / 24  hours) 

A ij  is  the  fraction  of  machine  type  i required  to  perform  each  operation  j, 

( = 1 / (PQij  * Y,ij)  if  Pij  is  equal  to  1) 

Clij  discounted  operating  cost  of  each  machine  of  type  i to  process  operation  j 

at  period  t, 

F 

" discounted  capitalized  procurement  cost  of  an  additional  machine  of  type 

i at  period  t, 
h 

tJ  discounted  cost  of  holding  one  unit  of  inventory  of  wafer  type  j at 
period  t. 
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Before  we  present  the  problem  formulation  we  define  the  following  production, 
inventory,  and  machine  requirement  variables: 

Xtij  is  the  number  of  units  of  wafer  type  j assigned  to  machine  group  i at 
period  t, 

I,i  is  the  number  of  units  of  wafer  type  j held  in  inventory  at  period  t, 

Qti  is  the  number  of  additional  machine  of  type  i needed  at  period  t. 


A comprehensive  mathematical  formulation  of  the  multi-period  alternative 
machine  type  capacity  allocation  problem  with  duplicated  machines  may  be  expressed  as 
follows: 


Problem  MPCAP 


EEE  c«/>®*«  + EE  v« 

te  T ie  I je  J teT  je  J 

+ EEF»e« 

re  T ie  P 

(1) 

I (t- 1 )j  X tij  ~ I tj  ~ V tj  ’ 

ie  I 

Vfe  T,  V;e  J 

(2) 

'ZP,„X,IJ<N„U,I, 

jeJ 

Vfe  T,  Vie  S 

(3) 

jeJ  t = 1 

Vfe  T,  Vie  P 

(4) 

Qti  > 0 and  integer , 

Vf  e r,  Vze  I 

(5) 

0 <Xtij<MPij, 

Vfe  T,  Vie  /,V/e  J 

(6) 
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For  a given  Uti  we  define  btij  = ptij  / Uti  and  replace  constraints  (3)  and  (4)  with: 


\/t  e T,  Vi  e 5 (7) 


t 


'Lb^Xq  <n„  + xe 


Vf  e T,  Vi  e P (8) 


7 = 1 


The  objective  function  of  the  formulation  is  to  minimize  operating  cost  of 
machines,  procurement  cost  of  additional  machines,  and  inventory  holding  cost. 
Constraints  (2)  ensure  that  the  demands  at  each  period  are  satisfied.  Constraints  (7) 
express  the  capacity  limits  on  the  "secondary"  machine  groups  and  constraints  (8)  express 
the  capacity  limits  on  the  "primary"  machine  groups  where  capacity  may  be  increased  by 
adding  new  machines.  We  prefer  to  express  machine  capacity  restrictions  under  two  sets 
of  constraints  to  emphasize  two  different  classes  of  machine  groups.  These  two 
constraints  may  as  well  be  formulated  as  one  constraint  set  similar  to  (8).  Constraints  (5) 
impose  integrality  and  non-negativity  restrictions  on  the  additional  machine  variables. 
Constraints  (6)  express  the  usual  non-negativity  conditions  on  assignment  variables  and 
ensure  that  operations  are  assigned  to  machine  groups  that  are  capable  of  processing 
them.  M is  a sufficiently  large  number. 

In  the  next  section,  we  discuss  a Lagrangian  relaxation  based  solution  procedure 
to  obtain  a good  lower  bound  on  the  minimum  value  of  problem  MPCAP  and  a feasible 


upper  bound. 
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3.3.  Description  of  the  Lagrangian-based  Heuristic  Solution  Method 

Lagrangian  relaxation  has  been  successfully  employed  by  numerous  researchers  in 
the  past  in  many  integer/mixed  integer  programming  applications.  The  approach  is  based 
on  the  observation  that  many  difficult  integer  programming  problems  can  be  viewed  as 
they  are  composed  of  two  types  of  constraints:  “nice”  constraints  and  “complicating” 
constraints.  Lagrangian  relaxation  is  formed  by  multiplying  the  complicating  constraints 
with  corresponding  penalty  costs  (Lagrangian  multipliers)  and  including  them  into  the 
objective  function.  By  dualizing  those  constraints,  we  obtain  a problem  that  is  easy  to 
solve  and  whose  optimal  value  gives  a lower  bound  on  the  optimal  value  of  the 
minimization  problem.  Although  the  method  does  not  guarantee  that  a solution  will  be 
found  for  all  problems,  it  is  often  worth  trying  since  it  is  so  simple 


Step  1 Initialize  Lagrangian  multipliers. 

Step  2 Update  the  lower  bound  on  the  minimum  value  of  the  objective 
function  by  solving  relaxed  model. 

Step  3 Compute  the  machine  groups'  workloads  and  inventories  from  the 
current  solution  and  update  Lagrangian  multipliers  with  respect  to 
violations  of  the  constraints. 

Step  4 Update  the  upper  bound  on  the  minimum  value  of  the  objective 
function;  if  a stopping  criteria  is  met,  stop;  otherwise  go  to  Step2. 


Figure  3.2.  Outline  of  the  heuristic  procedure 
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compared  to  other  methods.  The  reader  is  referred  to  Geoffrion  (1974)  and  Fisher  (1981) 
for  a detailed  study  on  the  Lagrangian  relaxation  technique. 

Our  solution  method  is  based  on  the  Lagrangian  relaxation  of  the  capacity 
constraints  (7),  (8),  and  a set  of  cumulative  machine  requirements  constraints,  which  is 
discussed  in  the  next  section.  An  overview  of  the  overall  solution  procedure  for  problem 
MPCAP  is  presented  in  Figure  3.2.  Through  Lagrangian  relaxation  of  the  above 
mentioned  constraints  in  Step  2 we  obtain  two  subproblems,  one  with  assignment  and 
inventory  variables,  the  other  with  additional  machine  requirement  variables.  The  first 
problem  further  decomposes  into  several  subproblems  that  are  easily  solved  by 
inspection,  as  is  the  latter.  The  solutions  of  these  problems  are  used  to  compute  a lower 
bound  on  the  minimum  value  of  the  objective  function. 

In  Step  3,  the  assignment  scheme  and  inventory  levels  resulting  from  Step  1 are 
determined  and  the  resource  workloads  are  computed  using  equations  (2)  - (4).  The 
violations  in  the  constraints  are  employed  to  update  the  Lagrangian  multipliers  applying  a 
subgradient  optimization  technique.  Finally,  the  upper  bound  on  the  minimum  value  of 
the  objective  function  is  computed  in  Step  4.  The  violations  of  capacity  constraints  (3)  are 
eliminated  through  an  operation  shifting  procedure  and  then  a machine  reduction 
procedure  is  applied  to  equations  (4). 

3.3.1.  Computation  of  the  Lower  Bound 

Before  we  proceed  with  the  Lagrangian  relaxation  of  the  problem  MPCAP  we 
first  consider  adding  a set  of  valid  inequalities  as  lower  bounds  on  the  total  number  of 
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additional  machines  needed  to  satisfy  the  cumulative  demand  from  periods  1 through  t, 
where  t = \ , ...  ,r: 

V/eT  (9) 

ieP  r= 1 z=  1 


These  constraints  may  provide  better  bounds  in  the  relaxed  problem.  Kt  are  the  minimum 
cumulative  additional  machine  requirements  to  meet  the  demands  of  periods  1 through  t 
and  are  best  obtained  by  solving  the  following  problem: 


Problem  MPCMR 


Min  X X Q„ 


teT  ieP 


s-t.  Z2X,.  >Xr,,  VteT.VjeJ 


r=l iel 


r=l 


and  (5)  - (8) 


and  then  setting 


t T 


*,=E££e2,’  t=u-,T 


ieP  r=l  z=l 


It  is  as  difficult  to  solve  problem  MPCMR  as  problem  MPCAP.  Therefore  we  propose  the 
following  set  of  uncapacitated  problems  that  can  be  easily  solved  by  inspection  as  an 
alternative  procedure  to  obtain  K{. 
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Problem  CMR, 


Min  ZZZV^h/ 


r=l  iel  jeJ 


\/jeJ 


r=l  iel 


r=l 


\/ieI,VjeJ 


where 


r=l  ie/  jeJ 


The  solution  procedure  for  Problem  CMR(,  t=  1, . . r is  as  follows:  for  each  period  r= 
1, . . t assign  the  whole  volume  of  operation  j to  the  fastest  machine  group  i where  it 
consumes  least  amount  of  resources  regardless  of  the  machine  group’s  capacity  but 
subject  to  its  capability.  Identify  the  number  of  additional  machines  required  in  each 
machine  group  i.  Round  up  the  difference  between  the  total  number  of  machines  needed 
and  total  number  of  machines  currently  available. 


Proposition:  Constraints  (9)  and 


ieP  r=l 


are  equivalent. 


r 


ieP 

t>\ 


ieP 
t>  2 


ieP 

t=r 
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= I'0,i+2>-ue2, 

ieP  ieP 

f>l  t>  2 


+ 


-+Ze 


ti 


ieP 

t=r 


= Z t(‘-t  + VQ„ 

ieP  r=l 


Our  solution  method  is  based  on  the  Lagrangian  relaxation  of  the  capacity 
constraints  (7),  (8),  and  cumulative  machine  lower  bounding  constraints  (10).  The  overall 
solution  procedure  for  problem  MPCAP  is  outlined  in  Figure  3.2.  The  formulation  of  the 
relaxed  problem  is  as  follows: 


Problem  LR(A,/u) 

Min  Z Z Z + LL,  K,  + Z Z V-,  + Z Z F„Q« 

teT  iel  jeJ  teT  jeJ  teT  ieP 

- Z Z Z KQ«  -Z  Z Z ftO  - r + O0rt  + c (11) 

teT  ieP  r=l  teT  ieP  r= 1 

s.t.  (2),  (5),  and  (6) 
where  X,/i  > 0 and  C is  constant: 

c = -ZZv+Zv 

teT  iel  teT 


The  term  involving  Lagrangian  multiplier  Xfi  and  Qn  in  the  objective  function  can 


alternatively  be  expressed  as  follows: 
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Z X Z Le„ 

/eT-  is?  r=l 

= X [^1/01/  + ^2;  (01;  + 02;)  + — + ^rifeli  + 02;  + •••  + 0r;)] 

ieP 

= Z K^li  + ^2;  + - + X?l/  + (^2;  + - + ^ri  )Qli  + - + ^r;0ri  ] 

ieP 

= Z Z kQm  +Z Z ^02,  +-+ZZ  *„0rt 

;'eP  r=l  ieP  r=2  ;'eP  r=r 

= 2 1 1 KQ* 

teT  ieP  r=t 


Using  a similar  expression  for  the  term  with  the  Lagrangian  multiplier  ju  the  objective 
function  is  reformulated  as  follows: 

III  + V*  + 1 1 V» 

/eT1  16/  ye/  teT  jeJ 


-zz 

(ef  ieP 


-K^  + tr 

T = t 


t + U/O 


0;;  + C 


(12) 


Lagrangian  problem  LR(X,/u)  nicely  decomposes  into  two  subproblems,  one  involving 
the  X(ij  and  Ifj  variables  only  and  the  other  involving  Qn  variables  only.  First  problem 
further  decomposes  into  n uncapacitated  multi-period  single-item  production  planning 
problem  which  can  be  easily  solved  by  inspection.  A procedure  to  solve  this  problem  is 


as  follows: 


set  all  Xnj , Ifj  to  zero 
for  t = 1,  r do 
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i*  = argminjej  {cUJ  J3tij  + Atibtij  } 
set  Xtl  *j  — Vtj 
endfor 

for  t = r,  2 do 

for  r = (t-1),  1 do 

If  {(cu*j  Pti*j  + ki*bti*j  >cri*/  Pri*j  + Ki  *bTi*j  + hrJ)  & (Xn*j  > 0)}  for  break 
else  z—  t-1 
endfor 
if  (z>  0)  do 

for  z =z,  . . ,t-l  do  I Tj  = Itj  + X(i*j 
Xft  *j  — Xr[  *j  + Xu  *j 
Xfi*j  = 0 
Itj -Itj -Vtj 

endif 

endfor 


Second  problem  can  also  be  solved  by  inspection:  if  the  coefficient  of  the  integer 
variable  Qti  is  positive  then  Qa  is  set  to  zero.  Otherwise,  Qa  would  be  infinite.  We 
impose  an  upper  bound  for  Qti  to  avoid  an  unbounded  solution. 


68 


We  use  a subgradient  optimization  procedure  to  update  the  Lagrangian 
multipliers.  Given  initial  values  of  A0  and  ju,  we  generate  a sequence  of  Lagrangian 
multipliers  Ak  and  /u k through  the  addition  of  a direction  vector  which  is  multiplied  by  a 
step  size  0 where  0 is  a positive  scalar.  Hence,  the  updating  of  Lagrangian  multipliers 
of  the  constraint  sets  (7)  and  (8)  is  done  according  to  the  equations 


A kti=max 


0;2f“’  + ek 


( \ 

X - ff. 

\JeJ  ) 


Vt  e T,  Vi  e S 


r 

f 

, VI 

AT 

1 

II 

■ae  ’>85 

0;A*_1  + ek 

-Nti~ZQn 

r=1  A 

Vt  £ T,  Vi  £ P 


and  Lagrangian  multipliers  of  the  constraints  (10)  is  updated  according  to  the  equation 


fikt  = max  < 


0 ; /j,k  1 + 6 


K,-'Z'L(t-T  + \)QII 


ieP  r=l 


, VteT 


The  step  size  0k  is  updated  at  each  iteration  k using  the  following  equation: 
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The  step  size  Ok  depends  on  the  parameter  s\o<s‘  <0.025  ),  on  the  gap  between  the 
current  lower  bound  'LB(Xk'\  //"')  and  the  estimated  minimum  value  of  the  objective 
function  of  the  relaxed  problem,  which  is  approximated  by  the  upper  bound  UB  obtained 
by  applying  a heuristic  method,  and  on  the  Euclidean  norm  of  the  deviations  in  the 
relaxed  constraints  (7),  (8),  and  (10).  The  sequence  Sk  is  determined  by  setting  5°  - 0.025 
and  by  dividing  5k  by  1.5  whenever  LB(XkA , // kA)  does  not  increase  after  a fixed  number 
of  iterations.  We  terminate  this  procedure  when  one  of  the  following  stopping  criteria  is 
met: 

1 . An  iteration  number  limit, 

2.  Maximum  gap  between  the  lower  and  upper  bounds, 

3.  A limit  on  the  value  of  the  Euclidean  norm  of  the  deviations. 

The  Lagrangian  problem  LR(A,ju)  has  the  Integrality  Property,  that  is  the  optimal 
solution  does  not  change  if  we  drop  the  integrality  restriction  on  Qa.  Thus,  Lagrangian 
relaxation  cannot  do  better  than  LP  relaxation.  However,  we  find  promising  to  use 
Lagrangian  relaxation  since  the  problems  we  consider  are  of  very  large  scale  and 
Lagrangian  relaxation  can  provide  good  lower  bounds  substantially  faster  than  the 
standard  LP  relaxation. 

3.3.2.  Computation  of  the  Upper  Bound 

At  each  iteration  of  the  Lagrangian  relaxation  the  upper  bound  is  computed  by 
first  modifying  the  solution  to  obtain  feasibility  and  then  by  applying  a procedure  to 
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improve  the  bound.  Since  the  Lagrangian  relaxation  results  in  a solution  that  is  feasible 
with  respect  to  the  inventory  balance  constraints,  we  must  ensure  that  the  production 
quantities  do  not  lead  to  over  usage  of  the  secondary  machines.  With  respect  to  the 
resource  constraints,  a feasible  production  plan  can  easily  be  created  by  shifting 
operations  from  overloaded  secondary  machine  groups  to  capacity  flexible  primary 
machine  groups.  The  selection  of  which  operations  to  be  shifted  is  done  by  sorting  the 
machine  group-operation  pairs  at  each  period  in  non-increasing  order  of  the  product  of 
consumption  and  operating  cost  coefficients  and  then  choosing  the  candidate  machine 
group-operation  pair  following  that  sequence.  Note  here  that  since  another  secondary 
machine  group  may  as  well  be  under  utilized  in  the  original  plan  and  therefore  may  be 
included  in  the  candidate  list.  We  continue  performing  this  shifting  procedure  until  the 
machine  requirements  in  all  overloaded  secondary  machine  groups  at  all  periods  are  at 
their  respective  capacity  levels.  An  operation  quantity  may  be  partially  reassigned  if  it  is 
not  necessary  to  shift  the  whole  production  volume  to  bring  the  total  resource 
consumption  in  certain  machine  group  to  capacity.  The  procedure  is  as  follows: 

for  t = 1,  . . r sort  OperCostdj*Consumptiondj  in  non-increasing  order 
for  t = 1, r do 

for  i = 1,  . . mj  do 

while  ( MachReqdd  > AvailMachp ) do 

shift  X(j  from  MachGroupi  to  MachGroupi * following  the  sequence 
update  MachReqdd  and  MachReqdd  * 
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endwhile 

endfor 

endfor 

for  t — 1,  ....  r do 

for  i = mj,  . . .,  m do 

while  ( MachReqdti  > Avail  Mach  a + MachUBtj)  do 

shift  Xy  from  MachGroupf  to  MachGroupi*  at  period  t or  t-1  following 
the  sequence 

update  MachReqdn  and  MachReqdti * or  MachReqd(t-j )i* 
endwhile 
endfor 
endfor 

MachUB  is  the  upper  bound  on  the  number  of  additional  machines  imposed  to  avoid 
unbounded  solution  and  my  is  the  number  of  secondary  machine  groups. 

After  a feasible  production  plan  is  constructed,  in  the  second  stage  we  try  to 
improve  the  solution  by  means  of  operation  transfers  between  primary  machine  groups 
within  the  same  period  or  different  periods.  An  operation  can  be  transferred  from  a period 
to  an  earlier  period  only  in  order  to  not  lead  to  any  backorders.  The  procedure  is 
performed  in  a similar  fashion  using  the  above  mentioned  sequence  of  machine  group- 
operation  pairs  starting  with  the  current  period  and  then  proceeding  with  the  preceding 
periods  in  an  attempt  to  avoid  inventories. 
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3.4.  Experimental  Analysis 

The  heuristic  solution  procedure  is  coded  in  C programming  language.  A series  of 
computational  results  is  carried  out  on  a PC  with  300  MHz  Pentium  II  processor  using 
several  sets  of  randomly  generated  problems  of  differing  sizes.  The  performance  of  the 
procedure  is  evaluated  using  various  combinations  of  cost,  utilization,  and 
primary/secondary  machine  groups  structures.  Procurement  cost  coefficient  F tu  operating 
cost  coefficient  Cftj,  and  utilization  Un  are  drawn  from  uniform  distributions  using  the 
parameters  shown  in  Table  3.2.  The  operating  cost  is  obtained  as  a percentage  of  the 
procurement  cost,  i.e.  operating  cost  is  on  the  average  20%  of  the  procurement  cost  using 
t/(0.1,  0.3),  40%  of  the  procurement  cost  using  1/(0. 1,  0.7),  etc.  Fn  and  cnj  are  generated 
for  the  first  period  and  discounted  for  succeeding  periods.  Similarly,  Utij  are  generated 


Table  3.2.  Experimental  design 


Parameters 

Values  used 

Total 

Procurement  cost 

(10,  11),  (10,  14),  (10,  17),  (10,  20) 

4 

Operating  cost3 

(0.1,  0.3),  (0.1,  0.7),  (0.1,  1.1),  (0.1,  1.5),  (0.1,  1.9) 

5 

Utilization 

(0.65,  0,75),  (0.65,  0.85),  (0.65,  0.95) 

3 

Number  of  periods 

20,  40 

2 

Number  of  machine  groups/operations 

10/200  15/300  20/400  25/500  30/600 

Machine  group  structures6 

29 

Total  number  of  problems 

3480 

a Variable  cost  is  a fraction  of  the  procurement  cost 
b 5 for  10/200  and  15/300,  6 for  20/400  and  25/500,  and  7 for  30/600 
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for  the  first  period  and  then  decreased  using  a certain  constant  ratio  for  the  following 
periods  considering  decline  in  the  utilization  of  machines  over  time.  Production  volumes 
and  production  capacities  are  drawn  from  uniform  distributions  U{  1,400)  and 
U( 50,  1500),  respectively.  The  number  of  process  steps  repeated  throughout  the 
manufacturing  of  the  wafer  is  also  generated  form  a uniform  distribution  U{ 3,  10). 

We  tested  our  procedure  on  ten  problem  types  over  a wide  range  of  values  of 
number  of  machine  groups  and  operations.  Each  problem  type  is  solved  for  different 
primary/secondary  machine  group  structures  with  the  same  input  data.  The  smallest 
problem  we  consider  consists  of  44,340  constraints,  44,060  variables,  of  which  60  are 
integers  and  the  largest  problem  has  747,160  constraints,  744,960  variables,  of  which  960 
are  integers.  A summary  of  the  results  is  presented  in  Table  3.3.  The  gap  between  the 
upper  bound  and  lower  bound  is  computed  as  (UB-LB)AJB.  Although  the  worst  case 
performances  do  not  seem  promising  in  terms  of  absolute  values,  we  were  able  obtain 
fairly  good  solutions  relatively  rapidly  for  large  problems.  Note  here  that  it  took  more 
than  80  minutes  to  obtain  even  the  LP  optimal  solution  in  some  problem  instances  of  40 
periods,  500  operations,  and  25  machine  groups  that  we  tested  using  the  CPLEX® 
optimization  software.  The  average  CPU  time  of  all  360  problem  instances  of  this 
structure  using  our  solution  procedure  is  337  seconds  and  the  maximum  CPU  time  is  520 


seconds. 
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Figure  3.3.  Mean  gaps  over  operating  costs  for  different  levels  of  utilization 
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Figure  3.4.  Mean  gaps  over  procurement  costs  for  different  levels  of  utilization 


76 


Figure  3.5.  Mean  gaps  over  utilizations  for  different  procurement  and  operating  costs 


Figure  3.3  illustrates  the  change  in  the  mean  gaps  with  respect  to  operating  costs 
for  different  utilization  distributions.  As  the  mean  of  the  distribution  of  operating  cost 
increases,  the  mean  gap  decreases.  On  the  other  hand,  no  significant  variation  is  observed 
in  the  mean  gap  with  respect  to  the  procurement  costs  (Refer  to  Figure  3.4).  Figure  3.5 
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Table  3.3.  Summary  of  results  showing  the  gaps  between  upper  and  lower  bounds 


Machine  Groups/Operations 


10/200 

15/300 

20/400 

25/500 

30/600 

Best 

0.0538 

0.0320 

0.0217 

0.0193 

0.0271 

20  Periods 

Average 

0.0809 

0.0565 

0.0499 

0.0459 

0.0515 

Worst 

0.1163 

0.0903 

0.0975 

0.0935 

0.0960 

Best 

0.0173 

0.0194 

0.0326 

0.0258 

0.0176 

40  Periods 

Average 

0.0541 

0.0319 

0.0487 

0.0406 

0.0327 

Worst 

0.1045 

0.0573 

0.0958 

0.0707 

0.0734 

exhibits  the  change  in  the  mean  gaps  with  respect  to  the  distribution  of  utilization.  We 
can  observe  that  the  solution  procedure  performs  better  as  the  mean  utilization  increases. 
This  is  due  to  the  fact  that  more  additional  machines  are  needed  if  the  utilization  is 


smaller.  Detailed  results  for  differing  problem  structures  are  included  in  Appendix  B. 


CHAPTER  4 

INTEGRATING  PRINTED  CIRCUIT  BOARD  SCHEDULING  AND  COMPONENT 

GROUPING  IN  AN  OPENSHOP  MANUFACTURING  ENVIRONMENT 

4.1.  Introduction 

Printed  circuit  board  (PCB)  assembly  consists  of  a PCB  and  various  electrical  and 
electronic  components  that  are  placed  onto  the  board.  There  is  a tremendous  variation  in 
the  technology  of  circuit  boards  and  the  components.  The  interested  reader  is  referred  to 
Kear  (1987)  and  McGinnis  et  al.  (1992)  for  a detailed  description  of  printed  circuit 
assembly  technology. 

Modem  PCB  production  typically  uses  computer  controlled  insertion  machines  to 
automatically  assemble  electronic  components  onto  the  boards  in  a flow  shop  type  of 
production  line.  Production  control  of  this  assembly  process  deals  with  several  problems 
such  as  the  selection  of  the  component  types  to  allocate  to  each  machine,  the 
determination  of  the  sequence  of  PCBs  for  production,  the  evaluation  of  the  arrangement 
of  feeders  onto  the  machine,  and  the  determination  of  the  sequence  of  the  placement 
operations.  A comprehensive  review  of  the  literature  in  the  PCB  assembly  research  area 
may  be  found  in  McGinnis  et  al.  (1992)  and  Askin  et  al.  (1994). 

We  discuss  the  component  family  formation  and  scheduling  problem  in  three 
stages.  The  first  stage  deals  with  the  problem  of  component  grouping  in  order  to 
minimize  the  number  of  machines  visited.  The  second  stage  loads  component  families  to 
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automatic  placement  machines  and  attempts  to  balance  the  workload.  The  third  stage 
addresses  the  scheduling  of  the  PCBs  to  machines  in  order  to  minimize  the  maximum 
completion  time  (makespan)  with  the  secondary  objective  of  reducing  the  mean  flow 
time.  Minimizing  maximum  completion  time  would  either  increase  the  production 
capacity  or  reduce  work-in-process  inventory  and  reducing  the  mean  flow  time  would 
also  reduce  work-in-process  inventory  without  adversely  affecting  production  capacity. 
Thus,  the  issues  addressed  are 

1 . Grouping  component  types  into  families, 

2.  Loading  component  families  to  machines,  and 

3.  Routing  PCBs  to  machines. 

In  this  study,  an  assembly  system  consisting  of  identical  decoupled  insertion 
machines  is  examined.  The  machines  place  electronic  components  onto  PCBs  using 
surface  mount  technology.  Each  PCB  has  a bill  of  components  specifying  the  component 
types,  the  number  of  each  type  of  component,  and  the  feeder  slots  that  each  component 
requires  on  each  machine.  For  the  environment  studied,  insertion  time  is  assumed  to  be  a 
fixed  (unitary)  constant  for  each  component  regardless  of  PCB  type  and  the  component’s 
location  on  the  board.  Thus,  the  total  insertion  time  of  all  components  on  a PCB  the  sum 
of  the  total  number  of  components  it  requires. 

We  also  assume  zero  machine  setup  times  when  switching  from  one  PCB  type  to 
another.  This  assumption  is  based  on  the  fact  that  components  are  already  loaded  on  each 
machine.  To  solve  the  allocation  problem  of  the  component  types  to  machines, 
component  families  are  formed  explicitly  considering  the  machine  capacity  constraint  in 
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terms  of  feeder  slots.  Note  here  that  feeder  slots  that  each  component  type  requires  may 
be  different  and  each  component  is  assigned  to  a single  machine.  When  a machine  starts 
placing  a component  type  on  a board,  no  interruption  is  allowed  until  the  whole  operation 
is  completed  for  that  board  (i.e.,  no  preemption  occurs).  Since  there  are  no  precedence 
requirements  on  the  operations  performed  PCBs  can  visit  the  machines  in  any  order. 
Furthermore,  if  the  machines  are  decoupled,  then  shop  configuration  would  be  an 
openshop. 

The  remainder  of  this  chapter  is  organized  as  follows.  Section  4.2  reviews  the 
relevant  literature  to  PCB  manufacturing  and  presents  an  overview  of  the  openshop 
scheduling.  Section  4.3  describes  the  formulation  of  the  clustering  problem  under 
consideration  and  discusses  the  algorithms  developed  for  component  family  formation, 
component  loading  and  PCB  scheduling.  Section  4.4  contains  the  experimental  analysis 
of  the  proposed  methodology  on  industry  data. 

4,2,  Relevant  Literature 

Over  the  last  two  decades,  many  researchers  have  studied  the  assembly  of  printed 
circuit  boards  in  an  effort  to  develop  models  to  support  the  process  planning,  production 
planning,  and  scheduling  decisions.  Typical  objective  in  these  decision  models  includes 
minimizing  cycle  time  and/or  setup  time.  Drezner  and  Nof  (1984),  Ball  and  Magazine 
(1988),  Ahmadi,  Grotzinger,  and  Johnson  (1988),  Francis,  Flamacher,  Lee,  and  Yeralan 
(1989),  Grotzinger  (1992)  address  the  problem  of  assigning  feeders  to  machines  and 
sequencing  placement  activities  to  minimize  the  cycle  time  for  each  board  type.  These 
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papers  also  discuss  issues  such  as  the  requirement  to  change  nozzles  for  some  component 
types,  concurrent  movement  of  pick  and  place  activities,  and  multiple  delivery  heads. 

Lofgren  and  McGinnis  (1986),  Carmon,  Maimon,  and  Dar-El  (1989),  Maimon 
and  Shtub  (1991)  consider  clustering  techniques  for  components  and  boards  to  reduce  the 
setup  time.  Lofgren  and  McGinnis  (1986)  group  similar  component  types  and  then  assign 
these  groups  to  balance  the  workload.  Ammons,  Lofgren,  and  McGinnis  (1985)  address 
the  problem  of  allocating  components  to  workstations  to  balance  workloads  and  to 
minimize  workstation  visits  in  an  environment  where  an  assembly  visits  a workstation 
only  once.  Workload  balancing  has  also  been  studied  by  Fathi  and  Taheri  (1989)  and 
Rajan  and  Segal  (1989)  with  a secondary  objective  of  minimizing  the  number  of 
machines  to  be  visited.  Crama,  Kolen,  Oerlemans,  and  Spieksma  (1990)  describe  a 
hierarchical  decomposition  and  mathematical  programming  approach  for  a line  with 
several  placement  machines  which  are  devoted  to  the  assembly  of  a single  product  with 
restrictions  on  the  component  types  that  can  be  loaded. 

In  a recent  paper,  Balakrishnan  and  Vanderbeck  (1999)  consider  a partial  setup 
strategy  of  mounting  frequently  used  components  permanently  on  each  machine  and 
present  a model  to  assign  product  families  to  parallel  assembly  lines  incorporating  both 
workload  balancing  and  setup  time  minimization  objectives.  Partial  setup  strategies  have 
been  discussed  by  Leon  and  Peters  (1996)  in  the  case  of  a single  placement  machine  to 
minimize  the  makespan  and  by  Peters  and  Subramanian  (1996)  in  the  case  of  multiple 
placement  machines  operating  in  parallel  to  balance  the  tradeoff  between  processing  time 
and  changeover  time. 
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Ahmadi,  Grotzinger,  and  Johnson  (1988)  and  Tang  (1988)  propose  the  use  of 
duplicate  reels  of  each  component  type.  Additional  feeder  tapes  of  the  same  component 
can  be  placed  onto  different  machines  to  balance  the  workload  or  onto  the  same  machine 
to  improve  the  cycle  time.  Although  this  practice  may  improve  the  solution  in  certain 
cases,  it  requires  additional  investment  in  parts,  planning,  and  setup. 

Lofgren,  McGinnis,  and  Tovey  (1991)  discuss  the  component  allocation  problem 
with  precedence  constraints  and  prove  that  the  general  problem  of  minimizing  the  number 
of  workstation  visits  is  NP-hard.  They  also  show  that  standard  heuristics  have  arbitrarily 
bad  worst  case  performance.  Askin,  Dror,  and  Vakharia  (1994)  consider  the  case  where 
precedence  constraints  do  not  exist.  They  develop  a methodology  to  group  PCBs  into 
production  families,  allocate  component  types  to  machines,  divide  families  into  board 
groups  with  similar  processing  times,  and  schedule  the  groups  in  order  to  minimize  the 
makespan  and  reduce  the  mean  flow  time.  In  an  environment  where  precedence 
constraints  on  the  order  of  component  placements  do  not  exist,  the  machine  scheduling 
problem  can  be  studied  in  the  context  of  an  openshop  scheduling  problem.  Since  our 
problem  setting  is  similar,  we  overview  the  openshop  scheduling  literature  next. 

The  openshop  machine  scheduling  problem  is  defined  as  follows:  n jobs  (/,,  ...., 

Jn ) have  to  be  processed  by  m machines  (M„  ....  ,Mm  ).  Job  Jj,  i= 1,  ...,  n,  requires  at 
most  m operations  (O/,,  .... , Olfn),  where  Oir  is  the  operation  performed  on  machine  r for 
job  i.  Operation  0[r  has  to  be  processed  (uninterrupted  if  preemption  is  not  allowed)  for 
t[r  time  units  on  machine  Mr,  r= 1,  ...  , m.  tir  is  a nonnegative  integer.  In  an  openshop, 
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there  is  no  restriction  on  the  order  in  which  the  operations  of  a job  are  to  be  performed. 
Each  job  can  be  processed  on  at  most  one  machine  at  a time  and  each  machine  can 
process  at  most  one  operation  at  a time. 

In  this  section,  the  three-fold  notation  a\  /3\y  of  Lawler,  Lenstra,  Rinnooy  Kan, 
and  Shmoys  (1989)  is  used.  The  first  field  a specifies  the  machine  environment,  the 
second  field  (3  indicates  the  job  characteristics,  and  the  third  field  prefers  to  the 
optimality  criterion.  In  our  case,  the  problem  is  indicated  as  Om  \ t(r  = 1 | Cmax  where  O 
stands  for  openshop  discipline  and  m for  the  number  of  machines.  If  the  number  of 
machines  is  arbitrary  then  m is  omitted.  Q denotes  the  completion  time  of  job  7/  and 
Cmax  = maxi  Ci  is  the  completion  time  of  the  last  job  in  the  schedule. 

Dror  (1992)  examines  openshop  scheduling  problem  with  machine  dependent 
processing  times  and  provides  a polynomial  time  algorithm  for  O \ tir  = tr , n > m | 

Cmax  and  O \ t[r  = tr , tx  > 2t2  > .. . > 2 tm , '^=2tr  ^ H I C/ . He  also  proves  that 
O | tir  = tr , m > 3,  n < m \ Cmax  is  NP-hard  and  gives  a linear  time  algorithm  for  02  | 

Hr  = tr  I XQ  • Vakharia  and  Qatay  (1997)  show  that  Dror’s  (1992)  algorithm  is  not 
optimal  and  propose  a refinement  of  the  procedure  while  proving  its  optimality.  Askin, 
Dror,  and  Vakharia  (1994)  address  the  PCB  family  formation  and  scheduling  problem  in 
an  openshop  environment  with  a primary  objective  of  minimizing  the  makespan  for 
assembling  a batch  of  boards  and  a secondary  objective  of  reducing  the  mean  flow  time. 

In  our  study,  tj/s  (i.e.,  insertion  time  for  PCB  i on  machine  Mr)  are  estimated  by 
the  component  requirements  of  the  PCBs  based  on  the  assignments  of  the  component 
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families  on  machine  Mr.  Hence,  the  problem  is  that  of  assigning  component  types  to  the 
machines  (to  estimate  t[r)  and  scheduling  the  PCBs.  This  scheduling  environment  can  be 
characterized  in  the  context  of  a flexible  openshop  since  all  insertion  machines  are 
identical  and  have  the  flexibility  of  processing  any  PCB  type.  In  that  sense,  the  machines 
become  distinct  only  after  the  component  families  are  assigned.  Thus,  family  formation 
plays  a critical  role  in  this  setting  where  the  grouping  of  component  types  dictates  the 
type  of  schedule  generated.  Further,  the  completion  time  of  jobs  is  determined  by  the 
assignments  of  families  to  machines  since  the  processing  time  is  a function  of  the  number 


Objective 


Stage  I 


Stage  II 


Stage  III 


Maximize  similarity  between 
components  subject  to  maximum 
slot  capacity  for  each  family 


Load  component  families  on 
machines  balancing  the  workloads 


Minimize  makespan  while  reducing 
mean  flow  time 


Figure  4. 1 . Overview  of  the  proposed  method 
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of  components  to  be  loaded  on  each  board.  In  order  to  achieve  a short  makespan, 
workload  balancing  between  machines  needs  also  to  be  considered  when  assigning  the 
component  families.  A lower  bound  on  the  makespan  is  given  by  the  maximum  amount 
of  components  to  be  placed  by  any  machine  and  this  amount  is  minimized  in  the  case  of 
perfectly  balanced  machines. 

4.3.  An  Integrated  Method  for  Component  Family  Formation.  Component  Loading,  and 

PCB  Scheduling 

In  this  section,  we  present  a three  stage  method  to  group  components  into 
families,  to  load  component  families  on  insertion  machines,  and  to  schedule  PCBs.  An 
overview  of  this  method  is  shown  in  Figure  4. 1 . In  the  first  stage,  component  families  are 
created  in  order  to  maximize  the  similarity  between  components  without  exceeding 
maximum  slot  capacity  for  each  family.  In  the  second  stage,  we  load  component  families 
on  the  insertion  machines  attempting  to  balance  the  workloads.  In  the  third  stage,  we 
schedule  PCBs  to  machines  in  order  to  minimize  the  makespan  with  the  secondary 
objective  of  reducing  the  mean  flow  time.  We  now  proceed  to  describe  details  of  each 
stage. 

4.3.1.  Stage  I:  Component  Family  Formation 

A comprehensive  0-1  integer  programming  formulation  of  a component  family 


formation  problem  is  as  follows: 
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max  Y,lLSijxij 

iaR  /€  P 

(1) 

CO 

'■ — 

ii 

(2) 

= 1> 

jeP 

i e R 

(3) 

I ft  xu  * C, 

ieR 

JeP 

(4) 

VI 

ieR,  j^P 

(5) 

Xij.yj  e{0,l}, 

ieR,jeP 

(6) 

where  i,j  are  indices  for  component  types,  R and  P are  the  set  of  component  types  and 
medians,  respectively.  Sy  is  the  similarity  coefficient,/?  is  the  number  of  medians 
(number  of  families)  to  be  identified,^-  is  the  number  of  feeder  slots  required  by 
component  type  i,  and  C is  the  capacity  of  each  machine.  The  decision  variables  are 

J 1 , if  component  i is  assigned  to  family  j e P 

[ 0,  otherwise 

[ 1 , if  component  j is  the  median  for  family  j e P 
■)'/  [ 0,  otherwise 

The  objective  function  (1)  maximizes  the  sum  of  similarities  between  each  pair  of 
component  types.  Constraint  (2)  sets  the  number  of  component  families;  constraint  set  (3) 
ensures  that  all  component  types  are  assigned  to  a single  family;  constraint  set  (4) 
enforces  the  feeder  slot  capacity  on  the  machines  and  constraints  (5)  ensure  assignment 
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of  a component  to  a family  that  is  created.  Thus,  the  above  clustering  problem’s  objective 
is  to  construct  p families  that  are  as  homogeneous  as  possible,  while  not  violating  the 
capacity  requirements  of  the  machines. 

Since  the  optimal  0-1  integer  programming  formulation  described  above  is  not 
practically  tractable,  we  propose  an  efficient  grouping  algorithm  to  derive  good  solutions 
to  this  problem.  An  overview  of  the  three  major  steps  in  our  procedure  is  as  follows: 

• First,  a set  ofp  medians  is  obtained  using  an  initial  heuristic.  Two  different 
heuristics  are  used  for  the  initialization  process. 

• Second,  component  types  are  assigned  to  their  most  similar  median  without 
violating  the  feeder  slot  capacity.  After  all  components  are  assigned,  the 
similarities  between  each  component  and  all  other  components  assigned  to  the 
same  median  are  computed.  The  set  of  medians  is  updated  if  a higher  similarity 
level  is  achieved  by  a component  other  than  the  median. 

• Finally,  after  the  above  process  is  complete,  all  possible  interchanges  between 
pairs  of  component  types  assigned  to  two  different  families  are  considered.  If  an 
interchange  can  improve  the  objective  function  value  without  violating  the 
capacity  constraint,  the  two  component  types  are  swapped.  The  procedure  is 
repeated  until  no  further  feasible  exchange  is  possible. 

Before  describing  each  step  of  the  procedure,  we  first  define  the  following  additional 


notation: 
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F 

number  of  component  families  (may  or  may  not  be 

prespecified) 

NF 

number  of  component  families  at  each  iteration 

Pik 

bP  most  similar  median  to  component  type  i,  k=  1,2,3,... 

CurrLoad(k) 

current  feeder  slots  required  by  all  components  types  assigned 
to  median  k 

FSi 

feeder  slots  required  by  component  type  i 

OBJ 

objective  function  value 

Pfemp 

temporary  set  of  medians 

SMS(P) 

sum  of  the  maximum  similarities  between  all  pairs  of  i eR  and 
JeP 

X max  StJ  for  P c R 

ieR  J*p 

~ ' i*j 

- oo  for  P = 0 

ssik 

sum  of  similarities  between  component  i and  all  components 
where  i is  assigned  to  median  k 

= X!  sij  xik  xjk  > i e R,ke  P,i  *k 

j*R 

j*i 

4.3. 1 . 1 . Initialization  Algorithms 

Two  heuristics  are  used  to  form  the  initial  set  of  medians.  The  first  is  a greedy 
algorithm  where,  for  a given  set  of  medians,  the  component  type  that  gives  the  greatest 
immediate  increase  in  the  value  of  the  sum  of  the  similarities  is  chosen  as  the  next 
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component  to  be  included  in  the  median  set.  In  the  second  algorithm,  the  pair  of 
component  types  that  has  the  greatest  similarity  are  merged  into  a new  component  type, 
provided  that  the  number  of  feeder  slots  required  by  both  does  not  exceed  the  machine’s 
capacity.  If  the  total  number  of  medians  is  prespecified,  the  algorithms  will  stop  when 
this  number  is  reached.  Otherwise  they  will  stop  when  the  value  of  the  sum  of  similarities 
fails  to  increase. 

Initial  Algorithm  1 (IA  I): 

StepO:  Let  P°=  0,  NF  = 0,  / =1 . 

Step  1:  Let  jt  = arg  max  jmKphl  SMS(Pt  l u {/}) 

Step  2:  If  SMS(PtA  u {/'})  < SMS(PtA),  stop.  PtA  is  an  initial  median  set. 

Otherwise,  set  Pt  = P1'  u {jj}  and  let  NF  = NF  + 1 . 

Step  3:  If  P{  = R,  or  if  NF  = F in  case  F is  prespecified,  stop.  R is  an  initial 

median  set.  Otherwise,  let  t = t +1,  go  to  Step  1. 

Initial  Algorithm  2 (IA  II): 

Step  0:  Let  P°  = 0 , NF  = |fl|,  t =1 . 

Step  I:  Identify  if  and  jt  such  that  Sy=  ma  x/<y  {Sy} 

Step  2:  If  Sy  = 0,  stop.  PtA  is  an  initial  median  set. 

Otherwise,  compute  total  number  of  feeders  required  by  it  and  jf  TF 
If  TF>  C let  Sy  = 0,  go  to  Step  5. 

Otherwise,  go  to  Step  4. 


Step  3: 
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Step  4:  If  NF  < F,  stop.  P1'  is  an  initial  median  set. 

Otherwise  go  to  Step  5. 

Step  5:  Set  NF  = NF-l;  merge  if  and  jf  to  create  a new  temporary  component  type 

kf.  Let  Pt  = P O u {jf} . Let  t = t +1 , go  to  Step  1 . 

4. 3.1.2.  Grouping  Algorithm 

After  initial  set  of  medians  is  determined  component  types  are  assigned  to  their 
most  similar  median  as  long  as  the  machine  slot  capacity  is  not  exceeded.  The  assignment 
procedure  is  as  follows:  for  each  component  type,  the  absolute  values  of  the  differences 
between  the  similarity  coefficients  of  the  component  and  its  first  and  second  most  similar 
medians  are  computed.  Then,  component  types  are  assigned  to  medians  in  non-increasing 
order  of  these  absolute  values.  By  doing  this,  the  assignment  of  a component  to  a very 
dissimilar  median  is  prevented. 

When  all  of  the  assignments  are  completed  for  each  component  type  in  every 
family  the  sum  of  similarities  between  that  component  and  other  components  assigned  to 
the  same  median  is  computed.  If  a component  type  other  than  the  median  achieves  the 
maximum  sum  of  similarities  then  the  set  of  medians  is  updated  by  adding  that 
component  type  which  improves  the  objective  function  value  and  removing  the  old 
median.  The  method  is  repeated  with  the  new  set  of  p medians  as  long  as  the  objective 
function  value  continues  to  improve  with  a new  component  type. 

The  details  of  the  procedure  to  group  components  into  families  are  as  follows: 
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Step  0:  Let  P'  = Pl  from  the  initial  algorithm  and  set  t =1. 

Step  1:  For  each  component  type  i,  sort  median / s in  non-increasing  order  of  Sy, 

i eR,j  eP. 

Step  2:  Compute  DIF[  = | Sipn  - SiPn  |.  Sort  DIF{ s in  non-increasing  order. 

Set  k = 1,  go  to  next  step. 

Step  3:  If  component  type  i is  not  assigned  and  CurrLoad(pik)  + FSj  < C, 

assign  i to  median  pip.  Let  CurrLoadip^)  = CurrLoadipip ) + FSi,  go  to 
Step  4.  Otherwise,  let  k = k +1 , go  to  Step  3. 

Step  4:  Set  Pjemp~  0 • Compute  SSjk’s. 

Step  5:  Find  the  greatest  SS^.  If  i is  not  a median,  replace  median  k with 

component  type  i:  i becomes  the  new  median.  Set  Pjemp  ~ P Temp  u {*}■ 
Go  to  Step  6. 

Step  6:  Compute  OBJ.  If  Ppemp  = P*-,  stoP- 

Otherwise,  let  Pl  = Pjemp , t ~ t +1,  go  to  Step  1 . 

4.3.I.3.  Improvement  Algorithm 

We  consider  all  possible  component  exchanges  between  two  families  to  further 
improve  the  similarity  between  components.  The  algorithm  is  as  follows: 


Step  0:  Let  ij  eR  and  k,l  eP 
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Step  1:  If  (Sik  + Sji ) > (* Sn  + Sjk)  and  CurrLoad(k ) + FS[  - FSj  < C and 

CurrLoad([)  + FSj  - FSi  < C,  switch  component  types  i and  j between 
median  k and  median  /. 

Step  2:  Update  CurrLoad(k ) and  CurrLoad(l). 

If  all  components  and  medians  are  considered,  stop. 

Otherwise,  go  to  Step  0. 

4.3.2.  Stage  II:  Component  Family  Loading 

Component  families  are  allocated  to  machines  attempting  to  keep  the  workload  on 
the  most  utilized  machine  as  low  as  possible.  First,  we  compute  the  processing  time 
needed  to  place  all  the  components  in  each  family  onto  the  PCBs  that  require  them.  Then 
component  families  are  loaded  on  machines  such  that  total  number  of  feeders  required 
does  not  exceed  slot  capacity  C and  maximum  of  the  total  processing  times  is  minimized 
to  achieve  a balanced  workload.  The  algorithm  is  as  follows: 

Step  0:  Compute  the  total  processing  time  of  each  component  family. 

Step  1:  Load  families  for  which  total  feeder  slot  requirements  are  at  capacity  C. 

(dedicating  the  machine  to  a unique  component  family) 

Step  2:  Set  the  upper  bound  on  the  makespan  to  the  largest  processing  time  among 


all  families. 
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Step  3:  Load  sequentially  the  remaining  component  families  to  the  next  available 

machine  up  to  the  upper  bound  on  the  makespan.  Skip  families  for  which 
the  machine  does  not  have  sufficient  feeder  slot  capacity. 

Step  4:  If  the  machine  is  loaded  up  to  its  slot  capacity,  repeat  Step  3 for  the 

remaining  machines. 

Step  5:  If  there  exist  unloaded  component  families,  increase  the  upper  bound  on 

the  makespan  by  the  smallest  processing  time  among  all  remaining 
families.  Go  to  Step  3. 

4.3.3.  Stage  III:  Scheduling  PCBs  on  Machines: 

Given  that  the  primary  objective  is  to  minimize  the  maximum  completion  time 
with  a secondary  objective  of  reducing  the  mean  flow  time,  PCB  scheduling  is  carried  out 
as  follows.  First,  we  assign  PCBs  to  machines  based  on  their  component  requirements. 
After  all  of  the  assignments  are  completed,  we  create  a PCB-Machine  Processing  Time 
matrix  where  each  entry  represents  the  time  based  on  the  number  of  components  each 
PCB  requires.  Next,  the  machines  are  sorted  in  the  non-increasing  order  of  their  total 

workload:  £ tn  ^ til  - • • - X Urn  • The  logic  behind  this  Longest 

Processing  Time  (LPT)  type  ordering  is  that  the  machine  currently  having  maximum 
workload  is  the  bottleneck  in  the  system  and  boards  are  scheduled  to  that  machine  first  so 
that  it  will  have  no  idle  time  (in  the  case  of  the  first  machine)  or  the  least  amount  of  idle 
time  (in  the  case  of  the  following  machines).  This  LPT  type  ordering  is  employed  in  an 
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effort  to  minimize  the  makespan.  Then,  the  boards  are  sorted  in  the  non-decreasing  order 
of  their  machine  processing  times.  This  Shortest  Processing  Time  (SPT)  type  ordering  is 
applied  to  reduce  the  mean  flow  time  and  PCBs  are  then  scheduled  on  each  machine.  If  a 
PCB  is  currently  assigned  to  a previous  machine,  it  is  scheduled  at  an  earlier  time  if  the 
machine  has  sufficient  idle  time.  Otherwise,  it  is  delayed  until  its  previous  operation  has 
been  completed. 

The  scheduling  algorithm  is  performed  as  follows: 

Step  0:  Assign  PCBs  to  machines  where  the  components  they  require  have  been 

loaded  and  create  PCB-Machine  Processing  Time  matrix. 

Step  1:  Sort  the  machines  in  non-increasing  order  of  the  total  processing  times. 

(LPT  rule  for  machines) 

Step  2:  Sort  PCBs  on  each  machine  in  non-decreasing  order  of  their  processing 

time  on  that  machine.  (SPT  rule  for  PCBs) 

Step  3:  Starting  with  the  machine  that  has  longest  processing  time,  i.e.  starting 

with  the  first  machine  in  the  sequence,  schedule  PCBs  in  the  SPT  order.  If 
a PCB  is  being  processed  on  a previous  machine,  schedule  it  to  an  earlier 
time  if  the  machine  has  sufficient  idle  time.  Otherwise,  delay  it  until  its 
previous  operation  has  been  completed. 
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4.4.  Experimental  Analysis 

4.4.1.  Industrial  Setting 

The  methodology  described  in  the  previous  section  to  form  component  families 
and  to  schedule  PCBs  to  the  insertion  machines  was  coded  in  Pascal  programming 
language.  The  data  were  obtained  from  a major  electronics  company  and  consist  of 
information  regarding  22  different  PCB  types  they  manufacture  and  290  distinct 
component  types.  The  data  include  component  types  required  by  each  PCB  type  along 
with  the  number  of  units  of  component  types  to  be  mounted  on  the  PCB  and  the  number 
of  feeder  slots  required  on  the  insertion  machines  to  hold  the  component  type.  The 
number  of  slots  is  3,  6,  or  9.  There  are  eight  automatic  surface  mount  placement  machines 
that  are  identical  and  the  total  number  of  feeder  slots  available  on  each  machine  is  1 80. 

We  first  reduce  the  component  set  by  combining  components  required  by  the 
same  set  of  PCBs.  This  helps  to  reduce  the  memory  requirements  and  run  time,  and  to 
prevent  unnecessary  data  handling.  Components  which  are  to  be  placed  on  the  same  PCB 
or  PCBs  are  merged  into  a new  component  type,  provided  that  the  sum  of  the  number  of 
slots  required  by  the  new  component  does  not  exceed  the  machine  capacity  of  180  slots 
that  are  available.  By  doing  so,  the  component  set  is  decreased  to  74  distinct  component 
types  and  slot  requirements  of  the  new  component  types  and  the  number  of  units  mounted 
on  each  PCB  are  updated.  Table  4.1  shows  the  information  for  PCB  type  #22  that 
originally  required  1 1 distinct  component  types  and  now  requires  1 0 after  the  reduction 


procedure. 
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For  each  i and  j,  i,j  e R,  and  i * j,  similarity  coefficients  (Sjj)  could  be  defined  in 
a variety  of  ways.  In  our  case,  we  defined  them  as  follows: 


Sf1" 


NCB, 


TNB.+  TNBj  - NCB0 


NCBij  is  the  number  of  PCBs  requiring  both  component  types  i and  j,  TNBj  is  the 
number  of  PCBs  requiring  component  type  i and  similarly  TNBj  is  the  number  of  PCBs 
requiring  component  type  j. 

beB 


UMib  is  the  number  of  units  of  component  i mounted  on  PCB  b,  UMjb  is  the  number  of 
units  of  component  j mounted  on  PCB  b,  and  B is  the  set  of  PCBs  requiring  either 
component  i or  j or  both. 1 


min{  UMib , UM  jb  } 
max{  UMib , UMjb } 
~NCB~ 


The  primary  reason  for  defining  similarity  coefficients  in  several  ways  was  to 
investigate  whether  the  families  identified  are  significantly  impacted  by  this  parameter 


1 Sfj  is  in  fact  a dissimilarity  coefficient.  In  that  case,  the  problem  becomes  a 
minimization  problem. 
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Table  4. 1 . Component  information  for  PCB  type  # 22 


Component  No. 

Units  Mounted 

Feeder  Slots  Required 

1 

1 

6 

3 

2 

6 

5 

1 

6 

13 

4 

9 

15 

3 

6 

27 

46 

6 

46 

20 

6 

60 

7 

6 

71 

1 

9 

74 

74 

12 

choice.  In  what  follows  a numerical  example,  illustrating  the  computation  of  the  three 
similarity  coefficients  discussed  earlier,  is  given.  Table  4.2  lists  PCBs  that  require 
component  46  and  component  60  along  with  the  number  of  component  units  to  be 
mounted.  Noting  that  the  number  of  PCBs  that  require  component  type  46  is  8,  the 


Table  4.2.  Data  on  PCB  types  requiring  components  #46  and  #60 


Component  No.  46 

Component  No.  60 

PCB  Tvoe 

Number  of  units 
to  be  mounted 

PCB  Tvne 

Number  of  units 
to  be  mounted 

8 

8 

11 

5 

10 

10 

14 

8 

11 

3 

15 

18 

14 

2 

16 

8 

16 

2 

18 

1 

20 

14 

19 

1 

21 

20 

20 

2 

22 

20 

21 

7 

22 

7 
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number  of  PCBs  that  require  component  type  60  is  9 and  the  number  of  PCBs  that  require 
both  component  types  is  6,  the  similarity  coefficients  between  these  two  component  types 
are  calculated  as  follows: 


J 46,60 


8 + 9-6 


= 0.545, 


0 0 3 2 0 2 

— h H 1 1 1 

8 10  5 8 18  8 


J 46,60 


0 2_  1_  1_ 

— + 14  + 20  + 20  =ft324 


\36.60  = (8  - 0)2  +(10-0)2  +(3-5)2  +(2-8)2  +(0-18)2  +(2-8)2  + 

(0-1)2  +(0-l)2  +(14-2)2  +{20 -l)1  +{20 -If  =1048 


4.4.2.  Grouping  Component  Types  and  Scheduling  PCBs 

In  order  to  identify  the  median  components  for  each  component  family,  the  two 
initialization  algorithms  (IA I and  IA II)  are  used.  The  algorithms  require  the  user  to 
prespecify  the  number  of  medians,  hence  the  number  of  families  to  be  formed.  In  the  first 
algorithm,  if  the  user  does  not  to  prespecify  number  of  component  families  the  algorithm 
stops  identifying  new  medians  when  the  sum  of  similarities  between  the  component  pairs 
in  each  family  reaches  to  maximum.  At  the  end,  both  algorithms  identify  a number  of 
components  that  will  be  used  as  the  initial  median  set  in  the  grouping  algorithm.  Random 
generation  of  components  as  initial  medians  was  also  considered.  However,  heuristics 
IA  I and  IA  II  provided  better  performance  both  in  terms  of  quality  of  the  results 
measured  by  the  similarity  criterion  and  execution  time. 
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After  the  initial  medians  are  identified,  grouping  and  improvement  algorithms 
determine  component  families  for  the  specified  number  of  medians  as  described  in  Stage 
I.  A large  range  of  families  (8  to  40)  has  been  examined  for  experimental  purposes.  The 
reason  for  doing  so  is  that  the  component  family  population  and  number  of  families  may 
have  a significant  affect  on  the  assignment  of  components  to  machines,  therefore  on  the 
schedule.  Since  the  processing  times  for  PCB  types  on  each  machine  differ  significantly 
and  component  family  formation  algorithms  do  not  require  any  intensive  computations  it 
may  be  useful  to  investigate  different  number  of  families  to  compare  the  performance 
results. 

After  all  component  families  are  identified,  we  compute  total  number  of 
component  units  to  be  mounted,  i.e.,  total  processing  time,  for  each  family.  Next,  all 
families  requiring  1 80  feeder  slots  are  loaded  to  insertion  machines.  Then,  we  find  the 
maximum  processing  time  among  all  component  families  as  an  upper  bound  on  the 
workload  on  each  machine  and  start  loading  remaining  families  in  the  order  of  their 
formation  up  to  that  upper  bound.  The  same  procedure  is  repeated  gradually  increasing 
the  upper  bound  by  the  smallest  processing  time  among  all  remaining  families.  Regarding 
the  capacity  of  machines  note  that  total  number  of  feeder  slots  required  by  all  components 
is  1208  compared  to  a total  capacity  of  1440  (8  machines  times  180)  slots. 

Finally,  when  all  component  families  are  loaded,  a lower  bound  on  the  makespan 
is  obtained  as  the  total  processing  time  on  the  machine  mounting  the  maximum  number 
of  component  units  among  all  machines.  PCBs  are  then  assigned  to  machines  based  on 
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their  component  requirements  and  scheduled  following  the  algorithm  in  Stage  III 
described  in  the  previous  section. 

4.4.3.  Computational  Results 

The  charts  in  Figure  4.2,  4.3,  and  4.4  show  the  completion  times  and  mean  flow 
times  obtained  for  different  number  of  component  families  using  three  different  similarity 
coefficients  (SC).  Shorter  completion  times  are  achieved  when  the  number  of  component 
families  varies  between  14  and  20,  except  in  the  case  of  SCI  where  shortest  completion 
time  is  achieved  forming  36  component  families  using  IA  II;  when  IA I is  used  makespan 
stays  the  same  for  number  of  families  greater  than  12.  The  utilization  of  SC2  and  SC3 
results  in  shorter  completion  times  for  almost  all  numbers  of  families  compared  to  SCI. 
Furthermore,  SC2  outperforms  SC3  in  terms  of  both  makespan  and  mean  flow  time. 
Regarding  the  performance  of  the  Scheduling  Algorithm,  out  of  33  different  numbers  of 


Table  4.3.  Summary  of  results  for  IA  I 


IA  I 

Number  of 
Families 

Makespan 

Gap 

MFT 

Number  of 
Visits 

SC  1 

18 

1033 

15.97% 

499 

77 

SC  2 

14 

893 

2.80% 

443 

96 

SC  3 

13 

992 

12.50% 

465 

103 
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SCI 


Component  families 


Figure  4.2.  Makespan  and  mean  flow  time  using  similarity  coefficient  1 
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Figure  4.3.  Makespan  and  mean  flow  time  using  similarity  coefficient  2 
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SC  3 


Component  families 


Figure  4.4.  Makespan  and  mean  flow  time  using  similarity  coefficient  3 
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families  examined,  the  lower  bound  on  the  makespan  is  achieved  in  43  instances  using 
SCI  (10  with  IA I and  33  with  IA II),  57  using  SC2  (28  with  IA  I and  29  with  IA  II),  and 
64  using  SC2  (32  with  IA  I and  32  with  IA  II).  The  deviation  from  the  lower  bound  is 
1.14%  on  the  average  and  12.85%  in  the  worst  case.  The  results  obtained  using  two 
initialization  algorithms  are  summarized  in  Table  4.3  and  Table  4.4.  All  of  the 
completion  times  are  at  their  lower  bounds  for  the  component  families  created.  The 
percent  difference  between  the  makespan  found  and  the  minimum  makespan  is  indicated 
in  the  Gap  column.  Optimal  makespan  (of  868)  is  found  by  solving  the  following  Min- 
Max  Generalized  Assignment  Problem: 


mm  z 

s.t.  X tixU  ~ z’ 

ieP 

j gM 

'Zxu=l> 

jeM 

i eP 

X fi  xy  - c’ 

izP 

jeM 

Xij  G {0,1}, 

ieP,jGM 

where  M is  the  set  of  machines 

Component  family  formation  reduces  the  number  of  machines  to  be  visited 
substantially.  In  all  of  the  instances,  fewer  number  of  visits  is  achieved  compared  to  1 15 
total  machine  visits  (5.23  machines/PCB  on  the  average)  in  the  optimal  solution.  In  the 
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Table  4.4.  Summary  of  results  for  IA  II 


IA  II 

Number  of 
Families 

Makespan 

Gap 

MFT 

Number  of 
Visits 

SC  1 

22 

1106 

21.52% 

498a 

84 

SC  2 

15 

883 

1.70% 

425a 

94 

SC  3 

36 

894 

2.91% 

480 

83 

indicates  that  the  a smallest  MFT  is  obtained  among  all  component  family 
numbers  examined. 


case  where  the  shortest  makespan  is  obtained,  the  average  number  of  visits  is  4.27 
whereas  the  average  is  3.77  in  the  best  case. 

In  sum,  the  results  indicate  that  component  grouping  using  SC2  with  either  IA  I or 
IA  II  provides  near-optimal  solutions  and  better  number  of  machine  visits  for  the 
environment  studied.  The  grouping  and  scheduling  algorithms  developed  may  be  applied 
in  other  openshop  environments  where  component  mounting  times  do  not  depend  on  the 
location  of  the  component  feeders  on  the  machine,  no  setup  is  needed  when  switching 
from  one  PCB  to  another,  and  machines  are  identical  in  terms  of  placement  times  and 
number  of  feeder  slots  required  for  each  component. 


CHAPTER  5 

CONCLUSIONS  AND  DIRECTIONS  FOR  FUTURE  RESEARCH 


In  this  research,  we  first  presented  a mathematical  model  for  assigning  operations 
to  machine  groups  and  determining  additional  machine  requirements.  We  have  proposed 
a Lagrangian  relaxation/decomposition  based  solution  procedure  given  the  intractability 
of  the  model.  Our  approach  is  shown  to  provide  near  optimal  solutions  in  short  CPU 
times.  Our  methodology  may  be  applied/adapted  for  other  similar  manufacturing 
applications  where  resource  allocation  and  capacity  expansion  decisions  are  to  be  made. 

This  first  model  is  extended  to  consider  multi-period  production  planning  by 
including  inventory  variables  and  inventory  balance  constraints  between  consecutive 
periods  to  guarantee  consistency  of  production  volumes.  In  that  context,  we  propose  a 
mathematical  model  for  alternative  machine  type  capacity  allocation  problem  with 
machine  duplication  in  semiconductor  manufacturing  and  a Lagrangian-based  heuristic 
solution  procedure  to  solve  this  problem.  Our  approach  is  shown  to  provide  fairly  good 
feasible  solutions  in  a relatively  short  amount  of  time  without  necessitating  the  use  of  any 
optimization  tools.  The  solution  may  as  well  be  used  as  an  initial  solution  to  obtain  an 
optimal  or  better  solution  employing  an  optimization  software. 

The  method  may  be  applied  to  production  and  capacity  planning  problems  of 
other  industries  where  resource  allocation  and  new  equipment  acquisition  decisions  are  to 
be  made.  Demand  for  new  equipment  arises  primarily  from  two  sources:  replacement  of 
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existing  aged  equipment,  and  additional  equipment  required  to  meet  growth  in  demand 
for  the  firm's  products  and  services.  The  model  discussed  here  can  also  be  extended  to 
consider  all  processes  in  the  manufacture  of  a computer  chip  by  including  constraints 
between  consecutive  process  steps  to  guarantee  consistency  of  production  volumes. 
Furthermore,  backorders  may  also  be  considered  in  that  context.  The  approach  could  also 
be  extended  to  investigate  replacement  of  capacity  as  well  as  expansion  and  disposal  to 
adapt  to  arbitrary  changes  in  demand. 

A limitation  of  the  model  is  that  it  assumes  deterministic  technological  changes, 
i.e.,  all  technologies  available  in  the  future  are  assumed  to  be  known  at  the  beginning. 
Other  limitations  include  assumptions  of  space  availability  in  the  facilities  for  additional 
equipment  and  accurate  forecasts  of  future  costs  and  demands. 

Finally,  we  present  an  integrated  methodology  for  grouping  components  into 
families,  loading  these  families  of  components  on  individual  insertion  machines,  and 
finally,  scheduling  PCBs  on  the  machines.  Our  integrated  approach  is  shown  to  provide 
near-optimal  solutions  in  a reasonable  period  of  time.  Further,  when  comparing  it  to  a 
procedure  where  we  simply  load  individual  components  to  machines  rather  than  loading 
component  families  on  machines,  we  show  that  the  creation  of  families  results  in  fewer 
visits  of  PCBs  to  individual  machines  (leading  to  a reduction  in  materials  handling  costs). 

We  also  show  that  using  our  approach,  the  user  can  determine  the  "best"  number 
of  families  that  should  be  formed  based  on  scheduling  measures  rather  than  the  structural 
measures  used  in  prior  work  on  cellular  manufacturing.  This,  in  our  opinion,  is  a much 
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more  effective  approach  for  determining  the  number  of  families  since  it  looks  at  how 
actual  shop  performance  measures. 

Extensions  of  this  study  could  focus  on  several  aspects.  First,  more  effective 
algorithms  could  possibly  be  developed  for  solving  the  component  grouping  problem 
addressed  in  Stage  1 of  our  method.  Second,  although  we  have  only  experimented  with 
industry  data  in  PCB  manufacturing,  the  approach  could  also  be  extended/adapted  for 
other  similar  manufacturing  applications  such  as  machine  tool  manufacturing  which  is 
also  characterized  by  high  component  variety  and  few  final  products. 


APPENDIX  A 

EXPERIMENTAL  STUDY  RESULTS  FOR  OPERATION/MACHINE  TYPE 
ASSIGNMENT  METHODOLOGY 


The  results  of  the  experimental  analysis  discussed  in  Chapter  2 are  illustrated 
here.  The  numbers  on  the  upper  left  comer  of  each  table  indicate  number  of  machines  and 
number  of  operations,  respectively.  Z*  is  the  optimal  value  achieved  at  each  run,  LB  is 
the  lower  bound  on  the  minimum  value  of  the  objective  function,  and  UB  is  the  upper 
bound  on  the  minimum  value  of  the  objective  function.  Thus,  the  percent  value  is  the  gap 
between  the  optimum  value  and  the  lower  bound,  and  the  upper  bound  and  the  optimum 
value,  respectively.  The  average  is  also  shown  at  the  bottom  of  each  set  of  five  runs.  The 
final  table  displays  overall  means  of  six  problem  types.  This  table  excludes  problem  type 
with  40  machine  groups  and  400  operations  since  the  optimal  solutions  are  not  available. 
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110 


(4,92) 

Variable 

Cost 

Machine  Cost 

( 10,  12  ) 

(10,  15) 

( 10,20) 

( 10,  40  ) 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)/LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

0.003% 

4.003% 

0.346% 

5.338% 

0.479% 

3.217% 

1.978% 

0.000% 

0.008% 

0.052% 

0.016% 

0.110% 

1.575% 

0.035% 

1.188% 

0.000% 

(0.1,  0.3) 

0.006% 

0.130% 

0.305% 

0.000% 

0.031% 

0.165% 

2.103% 

0.000% 

0.004% 

0.000% 

0.050% 

0.000% 

0.602% 

0.019% 

1.528% 

2.076% 

0.014% 

0.048% 

0.026% 

3.784% 

0.113% 

0.000% 

1.479% 

0.000% 

0.007% 

0.847% 

0.149% 

1.846% 

0.560% 

0.687% 

1.655% 

0.415% 

0.01 1% 

2.989% 

0.089% 

0.007% 

0.936% 

1.664% 

0.019% 

0.000% 

0.016% 

0.000% 

0.006% 

0.005% 

0.022% 

0.449% 

0.068% 

0.000% 

0.60 

(0.1,  0.9) 

0.058% 

0.000% 

0.058% 

0.070% 

0.041% 

0.000% 

0.016% 

0.363% 

0.022% 

0.138% 

0.047% 

0.000% 

0.002% 

0.000% 

0.052% 

0.000% 

0.112% 

3.023% 

0.003% 

0.058% 

0.043% 

0.000% 

0.116% 

0.017% 

0.044% 

1.230% 

0.041% 

0.028% 

0.209% 

0.423% 

0.054% 

0.076% 

0.018% 

0.000% 

0.035% 

0.032% 

0.098% 

0.000% 

0.077% 

0.916% 

0.023% 

0.106% 

0.016% 

0.000% 

0.058% 

0.000% 

0.034% 

0.078% 

(0.1,  1.9) 

0.293% 

1.219% 

0.061% 

0.000% 

0.250% 

0.112% 

0.01 1% 

0.116% 

0.015% 

2.259% 

0.012% 

0.161% 

0.122% 

0.066% 

0.084% 

4.311% 

0.055% 

1.734% 

0.085% 

0.000% 

0.048% 

0.164% 

0.002% 

0.117% 

0.081% 

1.064% 

0.042% 

0.039% 

0.115% 

0.068% 

0.042% 

1.108% 

0.000% 

0.086% 

0.467% 

0.125% 

0.007% 

0.042% 

0.135% 

0.000% 

0.016% 

0.000% 

1.293% 

0.108% 

0.009% 

0.085% 

1.318% 

0.000% 

(0.1,  0.3) 

0.032% 

0.161% 

0.020% 

0.039% 

0.003% 

0.147% 

0.542% 

0.000% 

0.015% 

0.024% 

0.005% 

0.000% 

0.025% 

0.252% 

0.020% 

0.705% 

0.012% 

0.000% 

0.176% 

0.091% 

0.010% 

0.000% 

0.155% 

0.000% 

0.015% 

0.054% 

0.392% 

0.073% 

0.011% 

0.105% 

0.434% 

0.141% 

o 

■C 

0.085% 

4.320% 

0.002% 

0.331% 

0.025% 

0.382% 

0.012% 

0.435% 

5 

0.054% 

0.006% 

0.010% 

0.081% 

0.003% 

3.905% 

4.447% 

0.087% 

0.75 

(0.1,  0.9) 

0.007% 

0.000% 

0.059% 

0.000% 

1.108% 

0.000% 

0.037% 

0.752% 

0.056% 

0.000% 

0.056% 

0.179% 

0.107% 

0.000% 

0.042% 

0.952% 

0.015% 

0.000% 

0.057% 

0.121% 

0.028% 

0.000% 

0.333% 

0.000% 

0.043% 

0.865% 

0.037% 

0.142% 

0.254% 

0.857% 

0.974% 

0.445% 

0.057% 

0.060% 

0.245% 

5.529% 

0.174% 

0.031% 

0.058% 

0.135% 

0.082% 

0.000% 

0.046% 

0.198% 

0.088% 

0.049% 

0.167% 

0.000% 

(0.1,  1.9) 

0.039% 

0.000% 

0.059% 

0.000% 

0.016% 

0.183% 

0.682% 

0.116% 

0.030% 

0.000% 

0.051% 

3.554% 

0.050% 

0.691% 

0.171% 

0.000% 

0.018% 

0.000% 

0.073% 

0.000% 

0.064% 

0.210% 

3.894% 

0.000% 

0.045% 

0.012% 

0.095% 

1.856% 

0.079% 

0.233% 

0.994% 

0.050% 

0.042% 

0.113% 

0.002% 

0.036% 

0.000% 

0.000% 

0.000% 

0.000% 

0.012% 

7.437% 

0.033% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

(0.1,  0.3) 

0.001% 

0.121% 

0.028% 

0.009% 

0.001% 

0.057% 

0.032% 

0.000% 

0.018% 

0.208% 

0.005% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

0.021% 

0.000% 

0.017% 

0.114% 

0.879% 

0.000% 

0.027% 

0.000% 

0.019% 

1.576% 

0.017% 

0.032% 

0.176% 

0.011% 

0.012% 

0.000% 

0.102% 

0.020% 

0.017% 

0.000% 

0.160% 

0.100% 

0.043% 

0.048% 

0.076% 

5.781% 

0.059% 

0.000% 

0.008% 

0.462% 

0.019% 

0.000% 

0.90 

(0.1,  0.9) 

0.011% 

0.000% 

0.034% 

0.000% 

0.108% 

0.000% 

0.024% 

0.000% 

0.010% 

0.174% 

0.020% 

0.296% 

0.090% 

0.402% 

0.032% 

0.000% 

0.015% 

0.000% 

0.095% 

0.210% 

0.042% 

0.000% 

0.000% 

0.000% 

0.043% 

1.195% 

0.045% 

0.101% 

0.082% 

0.193% 

0.023% 

0.010% 

0.064% 

0.000% 

0.067% 

0.000% 

0.017% 

0.000% 

0.042% 

0.796% 

0.008% 

0.000% 

0.209% 

3.888% 

0.051% 

0.035% 

0.239% 

3.548% 

(0.1,  1.9) 

0.098% 

4.132% 

0.003% 

0.043% 

0.136% 

0.085% 

0.182% 

0.000% 

0.060% 

0.152% 

0.050% 

4.028% 

0.164% 

0.016% 

0.052% 

0.000% 

0.091% 

0.039% 

0.135% 

0.000% 

0.009% 

0.472% 

0.049% 

0.610% 

0.064% 

0.865% 

0.093% 

1.592% 

0.076% 

0.122% 

0.113% 

0.991% 

Ill 


(4,120) 


0.60 


0.75 


0.90 


Machine  Cost 

Variable 

( 10,  12) 

( 10,  15  ) 

( 10,  20  ) 

( 10,  40  ) 

Cost 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

0.049% 

4.479% 

0.037% 

5.024% 

0.394% 

7.020% 

0.501% 

6.214% 

0.004% 

0.179% 

0.052% 

0.000% 

0.015% 

0.120% 

0.242% 

2.420% 

(0.1,  0.3) 

0.108% 

0.102% 

0.026% 

3.911% 

0.846% 

4.224% 

1.502% 

7.582% 

0.003% 

0.081% 

0.072% 

0.000% 

0.184% 

6.234% 

2.448% 

1.904% 

0.047% 

0.051% 

0.023% 

0.359% 

0.297% 

1.601% 

1 .484% 

1.170% 

0.042% 

0.978% 

0.042% 

1.859% 

0.347% 

3.840% 

1.235% 

3.858% 

0.016% 

0.000% 

0.016% 

0.000% 

0.014% 

0.000% 

0.759% 

0.492% 

0.112% 

0.050% 

0.135% 

1.903% 

0.047% 

3.973% 

0.246% 

2.008% 

(0.1,  0.9) 

0.066% 

0.294% 

0.102% 

1.507% 

0.119% 

2.053% 

0.173% 

1.895% 

0.035% 

0.172% 

0.098% 

0.000% 

0.221% 

0.000% 

1 .009% 

0.005% 

0.025% 

5.818% 

0.020% 

4.776% 

0.032% 

0.754% 

2.752% 

0.728% 

0.051% 

1.267% 

0.074% 

1.637% 

0.086% 

1.356% 

0.988% 

1.026% 

0.092% 

1.914% 

0.057% 

0.216% 

0.042% 

0.983% 

0.275% 

0.000% 

0.279% 

1.024% 

0.225% 

0.519% 

0.428% 

1.178% 

0.390% 

5.809% 

(0.1,  1.9) 

0.242% 

0.000% 

0.052% 

1.434% 

0.487% 

0.965% 

0.144% 

0.046% 

0.685% 

0.659% 

0.059% 

2.166% 

0.253% 

3.013% 

0.565% 

0.000% 

0.138% 

2.342% 

0.055% 

1.076% 

0.493% 

0.070% 

0.163% 

0.000% 

0.287% 

1.188% 

0.090% 

1.082% 

0.341% 

1.242% 

0.308% 

1.171% 

0.101% 

0.521% 

0.192% 

1.042% 

1.313% 

0.800% 

2.740% 

3.502% 

0.041% 

0.114% 

0.287% 

6.564% 

0.507% 

2.432% 

0.543% 

1.396% 

(0.1,  0.3) 

0.046% 

0.375% 

0.134% 

2.022% 

1.419% 

0.000% 

0.245% 

2.402% 

0.046% 

0.000% 

0.009% 

0.056% 

0.506% 

0.000% 

1.538% 

0.000% 

0.041% 

0.044% 

0.006% 

0.212% 

0.591% 

0.000% 

0.015% 

0.000% 

0.055% 

0.211% 

0.126% 

1.979% 

0.867% 

0.646% 

1.016% 

1.460% 

0.038% 

0.319% 

0.058% 

0.306% 

0.415% 

0.308% 

1.045% 

3.599% 

0.023% 

0.359% 

0.568% 

2.773% 

0.192% 

0.086% 

1.509% 

0.000% 

(0.1,  0.9) 

0.136% 

0.000% 

0.154% 

0.931% 

0.288% 

3.228% 

2.630% 

2.449% 

0.112% 

0.192% 

0.123% 

3.842% 

0.324% 

0.951% 

0.788% 

5.895% 

0.111% 

0.310% 

0.021% 

0.139% 

0.031% 

0.444% 

0.051% 

1.836% 

0.084% 

0.236% 

0.185% 

1.598% 

0.250% 

1 .003% 

1.204% 

2.756% 

0.062% 

0.000% 

0.001% 

0.000% 

0.186% 

2.252% 

0.990% 

2.282% 

0.037% 

0.371% 

0.012% 

0.215% 

0.215% 

3.088% 

0.080% 

0.716% 

(01,1.9) 

0.020% 

3.158% 

0.056% 

0.079% 

0.089% 

0.370% 

0.041% 

0.000% 

0.043% 

1.380% 

0.188% 

0.451% 

0.051% 

0.000% 

1.612% 

0.063% 

0.196% 

2.592% 

0.023% 

0.000% 

0.110% 

0.142% 

1.834% 

4.428% 

0.071% 

1.500% 

0.056% 

0.149% 

0.130% 

1.170% 

0.911% 

1.498% 

0.100% 

6.454% 

2.872% 

6.358% 

3.760% 

3.810% 

2.964% 

0.000% 

0.053% 

6.976% 

0.539% 

4.922% 

3.245% 

8.811% 

0.522% 

7.021% 

(0.1,  0.3) 

0.400% 

0.000% 

0.064% 

0.000% 

0.427% 

6.654% 

2.264% 

8.727% 

0.314% 

6.740% 

0.063% 

6.286% 

2.378% 

4.073% 

5.721% 

0.000% 

0.146% 

6.111% 

3.063% 

1.542% 

2.366% 

1.477% 

3.231% 

6.768% 

0.203% 

5.256% 

1.320% 

3.822% 

2.435% 

4.965% 

2.940% 

4.503% 

1.057% 

6.996% 

2.016% 

8.133% 

2.232% 

2.885% 

0.789% 

1 .966% 

1.197% 

0.861% 

1.138% 

2.288% 

2.046% 

2.527% 

0.930% 

0.015% 

(0.1,  0.9) 

1.363% 

1.141% 

2.087% 

1.256% 

2.681% 

1.032% 

1.309% 

1.005% 

1.184% 

0.400% 

0.732% 

6.372% 

1.963% 

0.771% 

1.128% 

1.730% 

0.635% 

2.298% 

2.658% 

6.259% 

1.635% 

1.169% 

1 .062% 

3.843% 

1.087% 

2.339% 

1.726% 

4.862% 

2.111% 

1.677% 

1 .044% 

1.712% 

1.830% 

3.599% 

1.443% 

3.655% 

1.152% 

0.000% 

1.168% 

0.808% 

0.995% 

3.670% 

1.833% 

0.514% 

1.437% 

0.000% 

1 .259% 

2.484% 

(0.1,  1.9) 

1.084% 

0.000% 

1.372% 

0.000% 

1.971% 

0.555% 

0.707% 

4.208% 

1.548% 

0.848% 

0.853% 

0.661% 

0.465% 

0.039% 

0.493% 

0.247% 

1.312% 

0.154% 

1.750% 

0.204% 

0.541% 

0.343% 

0.525% 

0.117% 

1.354% 

1.654% 

1.450% 

1.007% 

1.113% 

0.187% 

0.830% 

1.573% 
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(4,200) 


.2 

5 


Machine  Cost 

Variable 

( 10,  12  ) 

(10,15) 

( 10,  20) 

( 10,  40  ) 

Cost 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

0.027% 

6.347% 

2.167% 

6.545% 

0.067% 

5.096% 

1.678% 

5.384% 

0.007% 

0.525% 

0.011% 

1.214% 

0.015% 

0.472% 

0.004% 

0.000% 

(0.1, 0.3) 

0.052% 

0.678% 

0.005% 

0.276% 

0.611% 

5.510% 

0.035% 

0.576% 

0.011% 

0.552% 

0.005% 

0.028% 

0.000% 

0.000% 

0.004% 

0.000% 

0.002% 

0.057% 

0.124% 

4.744% 

0.028% 

0.000% 

0.000% 

0.000% 

0.020% 

1.632% 

0.462% 

2.561% 

0.144% 

2.216% 

0.344% 

1.192% 

0.027% 

0.134% 

0.130% 

3.579% 

0.013% 

0.000% 

0.018% 

0.227% 

0.011% 

0.042% 

0.078% 

0.115% 

1.123% 

5.895% 

0.036% 

3.745% 

(0.1,  0.9) 

0.031% 

0.323% 

0.015% 

0.012% 

0.033% 

0.256% 

1.861% 

0.000% 

0.000% 

0.192% 

0.007% 

0.018% 

0.001% 

0.235% 

0.031% 

0.000% 

0.026% 

0.069% 

0.992% 

2.240% 

0.040% 

0.151% 

0.001% 

0.000% 

0.019% 

0.152% 

0.244% 

1.193% 

0.242% 

1.307% 

0.389% 

0.794% 

0.032% 

0.038% 

0.040% 

0.203% 

0.039% 

2.973% 

0.005% 

0.139% 

0.028% 

0.092% 

0.007% 

0.181% 

0.012% 

0.410% 

1.427% 

0.205% 

(0.1,  1.9) 

0.014% 

0.104% 

0.066% 

0.000% 

0.010% 

0.000% 

0.261% 

1.909% 

0.003% 

0.000% 

0.041% 

0.022% 

0.011% 

0.000% 

0.050% 

0.000% 

0.042% 

0.153% 

0.014% 

2.832% 

0.845% 

1.309% 

0.030% 

0.196% 

0.024% 

0.077% 

0.034% 

0.648% 

0.183% 

0.938% 

0.355% 

0.490% 

0.013% 

0.159% 

0.000% 

0.000% 

0.002% 

0.449% 

0.009% 

0.000% 

0.024% 

0.347% 

0.001% 

0.131% 

0.000% 

0.000% 

0.000% 

0.000% 

(0.1,  0.3) 

0.005% 

0.000% 

0.017% 

8.447% 

0.000% 

0.000% 

0.001% 

0.091% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

0.015% 

0.683% 

0.000% 

0.000% 

4.968% 

0.000% 

0.000% 

0.000% 

0.011% 

0.238% 

0.004% 

1.716% 

0.994% 

0.090% 

0.002% 

0.018% 

0.019% 

0.036% 

0.166% 

5.912% 

0.069% 

5.877% 

0.008% 

0.000% 

0.116% 

6.416% 

0.487% 

5.129% 

0.500% 

5.313% 

0.000% 

0.000% 

(0.1,  0.9) 

0.024% 

0.200% 

0.258% 

5.557% 

0.011% 

0.264% 

0.000% 

0.000% 

0.018% 

0.000% 

0.015% 

0.000% 

0.006% 

0.000% 

0.043% 

4.432% 

0.348% 

5.320% 

0.016% 

0.000% 

3.812% 

0.000% 

0.050% 

0.000% 

0.105% 

2.394% 

0.188% 

3.320% 

0.880% 

2.291% 

0.020% 

0.886% 

0.010% 

0.000% 

0.961% 

2.747% 

0.014% 

0.000% 

0.005% 

0.000% 

0.158% 

4.370% 

2.887% 

0.462% 

0.808% 

3.227% 

0.768% 

0.000% 

(0.1,  1.9) 

0.010% 

0.000% 

0.078% 

4.018% 

0.035% 

3.900% 

0.912% 

0.662% 

0.103% 

4.402% 

0.021% 

0.000% 

0.146% 

3.929% 

0.164% 

3.102% 

1.035% 

2.848% 

3.208% 

0.359% 

2.599% 

0.000% 

0.012% 

0.155% 

0.263% 

2.324% 

1.431% 

1.517% 

0.720% 

2.211% 

0.372% 

0.784% 

2.974% 

5.726% 

1.640% 

9.538% 

0.000% 

0.000% 

4.583% 

3.682% 

0.122% 

0.266% 

0.007% 

1.112% 

5.749% 

0.000% 

9.898% 

0.000% 

(0.1, 0.3) 

0.003% 

0.000% 

0.000% 

0.000% 

3.954% 

1.045% 

6.406% 

0.000% 

0.007% 

0.108% 

4.436% 

6.533% 

5.851% 

2.073% 

0.000% 

0.000% 

2.255% 

8.619% 

0.412% 

5.627% 

7.585% 

0.436% 

10.603% 

0.000% 

1.072% 

2.944% 

1.299% 

4.562% 

4.628% 

0.711% 

6.298% 

0.736% 

2.645% 

7.203% 

5.183% 

1.998% 

0.641% 

8.078% 

4.236% 

2.703% 

3.157% 

4.345% 

0.290% 

7.796% 

0.000% 

0.000% 

4.154% 

0.000% 

(0.1,  0.9) 

0.424% 

8.885% 

0.000% 

0.000% 

6.571% 

0.000% 

2.962% 

0.000% 

0.358% 

8.219% 

6.977% 

0.000% 

5.251% 

0.000% 

0.000% 

0.000% 

0.781% 

7.328% 

1.225% 

6.991% 

0.015% 

0.000% 

5.687% 

0.000% 

1.473% 

7.196% 

2.735% 

3.357% 

2.496% 

1.616% 

3.408% 

0.541% 

2.688% 

2.272% 

1.483% 

4.909% 

1.993% 

3.551% 

4.547% 

0.406% 

3.495% 

8.225% 

0.806% 

5.210% 

0.026% 

5.567% 

0.172% 

0.000% 

(0.1,  1.9) 

4.205% 

1.574% 

0.215% 

4.973% 

4.625% 

0.000% 

2.006% 

0.000% 

0.499% 

4.531% 

4.998% 

0.118% 

0.253% 

4.777% 

0.000% 

0.000% 

1.451% 

3.998% 

4.257% 

0.000% 

4.441% 

0.000% 

0.508% 

0.000% 

2.468% 

4.120% 

2.352% 

3.042% 

2.267% 

2.779% 

1 .447% 

0.081% 
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(15,92) 

Variable 

Cost 

Machine  Cost 

( 10,  12 ) 

( 10,  15  ) 

( 10,  20 ) 

( 10,  40  ) 

(Z*-LB)/LB 

(UB-Z»)/Z* 

(Z*-LB)/LB 

(UB-Z*)  / Z* 

(Z*-LB)/LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

0.114% 

0.000% 

0.680% 

3.120% 

0.128% 

3.290% 

0.355% 

0.000% 

0.047% 

1.116% 

0.195% 

0.248% 

0.292% 

0.448% 

0.202% 

0.135% 

(0.1,  0.3) 

0.448% 

1.708% 

0.588% 

0.170% 

0.090% 

0.352% 

0.061% 

1.583% 

0.140% 

0.025% 

0.497% 

0.000% 

0.017% 

1.251% 

0.077% 

0.125% 

0.204% 

1.082% 

0.153% 

0.207% 

0.168% 

2.195% 

0.620% 

0.000% 

0.191% 

0.786% 

0.423% 

0.749% 

0.139% 

1.507% 

0.263% 

0.369% 

0.341% 

2.776% 

0.100% 

2.759% 

0.156% 

0.301% 

0.023% 

2.682% 

0.077% 

1.499% 

0.548% 

0.348% 

0.125% 

0.000% 

0.240% 

0.481% 

0.60 

(0.1,  0.9) 

0.233% 

0.952% 

0.172% 

5.110% 

0.178% 

0.669% 

0.106% 

0.568% 

0.214% 

1.486% 

0.437% 

0.742% 

0.154% 

0.982% 

0.211% 

0.000% 

0.174% 

3.908% 

0.189% 

1.237% 

0.202% 

0.402% 

0.179% 

2.954% 

0.208% 

2.124% 

0.289% 

2.039% 

0.163% 

0.471% 

0.152% 

1.337% 

0.276% 

2.817% 

0.240% 

2.405% 

0.217% 

2.013% 

0.207% 

0.000% 

0.173% 

2.328% 

0.189% 

0.853% 

0.237% 

2.472% 

0.078% 

0.806% 

(0.1,  1.9) 

0.268% 

3.434% 

0.089% 

0.873% 

0.167% 

4.320% 

0.198% 

0.590% 

0.176% 

1.632% 

0.284% 

0.337% 

0.147% 

3.531% 

0.286% 

2.265% 

0.179% 

0.710% 

0.302% 

0.213% 

0.155% 

2.827% 

0.136% 

3.370% 

0.215% 

2.184% 

0.221% 

0.936% 

0.185% 

3.033% 

0.181% 

1.406% 

0.145% 

0.000% 

0.128% 

1.200% 

1.315% 

0.029% 

0.117% 

3.178% 

0.971% 

1.357% 

0.236% 

0.263% 

0.776% 

0.754% 

1.152% 

0.000% 

(0.1,  0.3  ) 

0.677% 

0.000% 

1.202% 

0.525% 

1.320% 

0.000% 

1.261% 

2.114% 

0.381% 

4.069% 

1.376% 

0.396% 

1.304% 

7.857% 

0.301% 

0.374% 

1.359% 

1.284% 

0.916% 

0.474% 

1.233% 

1.478% 

1.095% 

0.659% 

0.706% 

1.342% 

0.772% 

0.571% 

1.190% 

2.024% 

0.785% 

1.265% 

O 

•-C 

0.154% 

2.367% 

0.139% 

3.188% 

0.044% 

2.353% 

0.252% 

0.474% 

1 

0.108% 

3.021% 

0.379% 

1.308% 

0.179% 

0.000% 

0.342% 

2.418% 

<- 

0.75 

(0.1,  0.9) 

0.228% 

8.465% 

0.338% 

1.082% 

0.576% 

2.700% 

0.246% 

0.000% 

0.417% 

0.234% 

0.538% 

1.087% 

0.576% 

1.039% 

0.193% 

0.000% 

0.134% 

0.000% 

0.872% 

0.835% 

1.002% 

1.197% 

0.061% 

1.554% 

0.208% 

2.817% 

0.453% 

1.500% 

0.475% 

1.458% 

0.219% 

0.889% 

0.308% 

0.000% 

0.315% 

0.986% 

0.111% 

0.440% 

0.548% 

1.417% 

0.555% 

4.917% 

1.350% 

2.695% 

0.238% 

7.229% 

0.341% 

0.276% 

(0.1,  1.9) 

0.992% 

0.156% 

0.548% 

1.951% 

0.966% 

0.000% 

0.872% 

1.774% 

0.145% 

0.465% 

0.253% 

2.648% 

0.120% 

0.161% 

0.002% 

0.000% 

0.354% 

0.398% 

0.234% 

1.679% 

0.308% 

0.568% 

0.036% 

0.000% 

0.471% 

1.187% 

0.540% 

1.992% 

0.349% 

1.680% 

0.360% 

0.693% 

0.025% 

0.348% 

0.259% 

0.000% 

0.303% 

2.625% 

0.000% 

0.000% 

0.078% 

0.074% 

0.000% 

0.000% 

0.167% 

0.549% 

1.214% 

2.628% 

(0.1,  0.3) 

0.551% 

5.912% 

1.291% 

1.114% 

0.042% 

0.000% 

0.239% 

1.885% 

0.459% 

4.932% 

0.729% 

2.849% 

1.139% 

0.513% 

2.153% 

1.317% 

0.602% 

5.307% 

0.148% 

0.000% 

0.000% 

0.000% 

0.287% 

0.000% 

0.343% 

3.314% 

0.485% 

0.792% 

0.330% 

0.737% 

0.779% 

1.166% 

0.054% 

1.431% 

1.710% 

5.082% 

2.602% 

0.000% 

0.369% 

0.047% 

1.323% 

4.017% 

1.391% 

0.658% 

0.482% 

3.252% 

0.516% 

0.000% 

0.90 

(0.1,  0.9) 

1.656% 

5.604% 

0.594% 

0.000% 

0.722% 

1.422% 

0.577% 

0.000% 

1.257% 

7.184% 

3.069% 

0.000% 

0.507% 

1.607% 

0.788% 

0.000% 

0.345% 

4.830% 

1.061% 

8.649% 

0.455% 

5.018% 

0.326% 

0.000% 

0.927% 

4.613% 

1.565% 

2.878% 

0.953% 

2.260% 

0.515% 

0.009% 

1.831% 

0.040% 

0.639% 

2.322% 

0.984% 

0.000% 

0.348% 

0.081% 

0.499% 

2.032% 

0.342% 

2.329% 

0.292% 

0.321% 

0.012% 

0.000% 

(0.1,  1.9) 

0.947% 

1.063% 

0.967% 

1.039% 

1.713% 

0.000% 

0.034% 

0.000% 

1.198% 

0.963% 

0.575% 

3.018% 

0.802% 

0.639% 

0.016% 

0.000% 

0.876% 

0.089% 

1.078% 

0.700% 

1.087% 

0.023% 

0.124% 

1.001% 

1.070% 

0.837% 

0.720% 

1.881% 

0.976% 

0.197% 

0.107% 

0.216% 
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(15,120) 


0.60 


0.75 


Machine  Cost 

Variable 

(10,12) 

(10,15) 

( 10,20) 

( 10,  40 ) 

Cost 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

0.123% 

0.000% 

0.193% 

3.126% 

0.069% 

5.135% 

0.062% 

0.349% 

0.054% 

3.246% 

0.346% 

0.336% 

0.178% 

0.274% 

0.045% 

0.000% 

(0.1,  0.3) 

0.069% 

2.747% 

0.188% 

2.136% 

0.123% 

1.081% 

0.717% 

0.019% 

0.315% 

0.969% 

0.076% 

0.000% 

0.106% 

1.854% 

0.254% 

0.209% 

0.024% 

3.178% 

0.024% 

1.490% 

0.340% 

2.807% 

0.203% 

0.000% 

0.117% 

2.028% 

0.165% 

1.418% 

0.163% 

2.230% 

0.256% 

0.115% 

0.180% 

1.237% 

0.337% 

0.906% 

0.071% 

0.777% 

0.085% 

0.029% 

0.110% 

2.708% 

0.188% 

0.232% 

0.155% 

0.925% 

0.341% 

0.402% 

© 

© 

0.216% 

2.003% 

0.219% 

0.253% 

0.101% 

0.016% 

0.067% 

0.932% 

0.098% 

4.622% 

0.299% 

2.374% 

0.146% 

0.504% 

0.072% 

1.249% 

0.033% 

1.111% 

0.114% 

1.254% 

0.104% 

0.663% 

0.153% 

0.300% 

0.127% 

2.336% 

0.232% 

1.004% 

0.115% 

0.577% 

0.143% 

0.582% 

0.330% 

0.234% 

0.124% 

0.624% 

0.513% 

2.527% 

0.052% 

0.000% 

0.078% 

1.626% 

0.221% 

1.251% 

0.091% 

2.487% 

0.093% 

0.231% 

(0.1,  1.9) 

0.225% 

0.410% 

0.115% 

2.868% 

0.111% 

1.765% 

0.175% 

0.957% 

0.161% 

0.529% 

0.194% 

1 .366% 

0.200% 

0.802% 

0.104% 

0.393% 

0.115% 

1.655% 

0.160% 

2.155% 

0.239% 

1.282% 

0.061% 

0.271% 

0.182% 

0.891% 

0.163% 

1.653% 

0.231% 

1.772% 

0.097% 

0.370% 

0.609% 

0.773% 

0.023% 

2.042% 

0.037% 

1.397% 

0.510% 

0.000% 

0.786% 

0.843% 

0.335% 

2.293% 

0.475% 

0.156% 

0.220% 

0.000% 

(0.1,  0.3) 

0.324% 

0.780% 

0.033% 

0.879% 

0.057% 

0.464% 

1.296% 

0.436% 

0.380% 

3.543% 

0.501% 

1.960% 

0.165% 

0.414% 

0.069% 

2.753% 

0.225% 

1.750% 

0.317% 

0.786% 

0.371% 

0.000% 

0.076% 

0.333% 

0.465% 

1 .538% 

0.242% 

1.592% 

0.221% 

0.486% 

0.434% 

0.705% 

0.138% 

1.730% 

0.260% 

3.834% 

0.204% 

1.478% 

0.074% 

0.000% 

0.199% 

2.462% 

0.215% 

1.263% 

0.373% 

1.288% 

0.562% 

0.000% 

(0.1,  0.9) 

0.199% 

2.373% 

0.202% 

0.438% 

0.134% 

2.256% 

0.142% 

2.401% 

0.715% 

0.786% 

0.203% 

1.809% 

0.193% 

3.967% 

0.391% 

0.134% 

0.209% 

0.380% 

0.195% 

0.493% 

0.605% 

3.240% 

0.266% 

0.000% 

0.292% 

1 .546% 

0.215% 

1.567% 

0.302% 

2.446% 

0.287% 

0.507% 

0.223% 

0.582% 

0.569% 

0.225% 

0.473% 

0.000% 

0.169% 

1.211% 

0.459% 

0.365% 

0.225% 

2.710% 

0.102% 

0.296% 

0.256% 

0.824% 

(0.1,  1.9) 

0.250% 

4.013% 

0.465% 

0.676% 

0.286% 

0.801% 

0.243% 

0.189% 

0.494% 

0.881% 

0.192% 

3.949% 

0.285% 

2.188% 

0.125% 

0.718% 

0.134% 

1.954% 

0.486% 

0.873% 

0.208% 

0.545% 

0.146% 

0.266% 

0.312% 

1.559% 

0.387% 

1.686% 

0.271% 

0.766% 

0. 1 88% 

0.641% 

0.498% 

4.238% 

0.479% 

0.647% 

0.823% 

0.439% 

2.850% 

0.044% 

0.001% 

0.000% 

0.126% 

4.675% 

1.712% 

2.590% 

1.031% 

1.988% 

(0.1,  0.3) 

2.413% 

0.023% 

0.003% 

0.000% 

1.390% 

0.000% 

0.599% 

1.699% 

0.482% 

1.229% 

0.449% 

4.190% 

0.001% 

4.844% 

0.314% 

0.494% 

1.815% 

7.135% 

0.008% 

0.000% 

2.614% 

1.973% 

1 .506% 

0.012% 

1.042% 

2.525% 

0.213% 

1.902% 

1.308% 

1.969% 

1 .260% 

0.848% 

0.554% 

0.000% 

0.842% 

0.656% 

0.346% 

2.822% 

0.457% 

0.058% 

0.305% 

3.742% 

0.510% 

2.531% 

2.303% 

0.572% 

0.134% 

0.068% 

(0.1,  0.9) 

0.336% 

4.086% 

0.090% 

4.223% 

1.529% 

0.000% 

0.065% 

2.204% 

0.443% 

4.407% 

1.772% 

0.883% 

0.300% 

0.481% 

0.115% 

0.000% 

0.161% 

0.283% 

2.523% 

1 .483% 

0.607% 

1 .540% 

0.754% 

0.163% 

0.360% 

2.504% 

1.147% 

1 .955% 

1.017% 

1.083% 

0.305% 

0.499% 

0.630% 

5.619% 

0.905% 

1.247% 

0.536% 

1.936% 

0.328% 

0.753% 

0.182% 

3.355% 

1.353% 

0.842% 

0.461% 

1.392% 

0.209% 

3.132% 

(0.1,  1.9) 

0.359% 

0.000% 

0.462% 

0.675% 

1.360% 

2.676% 

0.076% 

0.264% 

0.469% 

0.327% 

0.365% 

0.000% 

0.190% 

0.000% 

0.056% 

0.000% 

0.581% 

1.360% 

0.713% 

0.449% 

0.063% 

1.055% 

0.013% 

0.000% 

0.444% 

2.132% 

0.760% 

0.643% 

0.522% 

1.412% 

0.136% 

0.830% 
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0.60 


‘■a 

jj 


0.75 


0.90 


Machine  Cost 

Variable 

( 10,  12) 

( 10,  15  ) 

( 10,  20) 

( 10,  40 ) 

Cost 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

(Z*-LB)  / LB 

(UB-Z*)/Z* 

(Z*-LB)  / LB 

(UB-Z*)  / Z* 

0.026% 

2.124% 

0.079% 

0.000% 

0.184% 

0.570% 

0.424% 

0.000% 

0.057% 

0.134% 

0.077% 

0.586% 

0.036% 

1.715% 

0.073% 

0.024% 

(0.1,  0.3) 

0.386% 

0.789% 

0.160% 

0.443% 

0.156% 

0.000% 

0.121% 

0.467% 

0.049% 

0.103% 

0.306% 

2.219% 

0.030% 

0.539% 

1.227% 

0.037% 

0.366% 

1.061% 

0.078% 

0.277% 

0.041% 

0.374% 

0.000% 

0.000% 

0.177% 

0.842% 

0.140% 

0.705% 

0.089% 

0.640% 

0.369% 

0.106% 

0.111% 

0.181% 

0.106% 

0.282% 

0.140% 

0.475% 

0.198% 

0.151% 

0.049% 

4.914% 

0.078% 

1.034% 

0.094% 

1.166% 

0.129% 

0.217% 

(01,0.9) 

0.192% 

0.537% 

0.040% 

0.563% 

0.049% 

2.188% 

0.070% 

0.588% 

0.116% 

0.137% 

0.100% 

2.061% 

0.338% 

0.611% 

0.035% 

0.000% 

0.149% 

2.605% 

0.091% 

1.043% 

0.164% 

0.555% 

0.048% 

0.150% 

0.123% 

1.675% 

0.083% 

0.997% 

0.157% 

0.999% 

0.096% 

0.221% 

0.166% 

1.830% 

0.104% 

0.816% 

0.151% 

0.871% 

0.119% 

0.345% 

0.184% 

1.422% 

0.114% 

1 .967% 

0.207% 

0.334% 

0.112% 

0.542% 

(0.1,  1.9) 

0.162% 

1.914% 

0.107% 

0.175% 

0.172% 

0.408% 

0.142% 

0.724% 

0.095% 

1.186% 

0.098% 

1.298% 

0.096% 

0.494% 

0.135% 

0.273% 

0.262% 

1 .946% 

0.106% 

0.176% 

0.183% 

1.124% 

0.072% 

0.647% 

0.174% 

1.660% 

0.106% 

0.886% 

0.162% 

0.646% 

0.116% 

0.506% 

0.167% 

0.015% 

0.000% 

0.105% 

0.028% 

1.564% 

0. 1 80% 

0.000% 

0.489% 

2.274% 

0.353% 

2.378% 

0.059% 

5.441% 

0.000% 

0.000% 

(0.1,03) 

0.024% 

0.473% 

0.448% 

2.218% 

0.162% 

0.564% 

1 .342% 

0.686% 

0.005% 

3.433% 

0.178% 

0.405% 

0.004% 

0.035% 

1.486% 

0.147% 

0.100% 

0.000% 

0.396% 

2.070% 

0.214% 

0.000% 

0.652% 

0.000% 

0.157% 

1.239% 

0.275% 

1.435% 

0.094% 

1.521% 

0.732% 

0.167% 

0.127% 

1.673% 

0.282% 

3.128% 

0.035% 

1.718% 

0.267% 

0.887% 

0.167% 

2.715% 

0.087% 

1.755% 

0.227% 

0.510% 

0.384% 

0.729% 

(0.1,  0.9) 

0.374% 

0.330% 

0.096% 

0.672% 

0.036% 

2.051% 

0.104% 

0.112% 

0.062% 

2.539% 

0.182% 

2.003% 

0.063% 

0.088% 

0.753% 

0.235% 

0.621% 

0.403% 

0.141% 

0.769% 

0.248% 

1 .479% 

0.351% 

1.313% 

0.270% 

1.532% 

0.158% 

1.665% 

0.122% 

1.169% 

0.372% 

0.655% 

0.226% 

1.925% 

0.269% 

0.912% 

0.268% 

2.664% 

0.083% 

0.000% 

0.453% 

0.672% 

0.135% 

1.796% 

0.410% 

0.545% 

0.075% 

1.015% 

(0.1,  1.9) 

0.215% 

0.000% 

0.246% 

1.530% 

0.228% 

1.523% 

0.140% 

0.027% 

0.184% 

1.233% 

0.420% 

0.692% 

0.150% 

0.198% 

0.058% 

0.826% 

0.543% 

0.234% 

0.271% 

1.344% 

0.114% 

1.858% 

0.082% 

0.000% 

0.324% 

0.813% 

0.268% 

1.255% 

0.234% 

1.358% 

0.088% 

0.374% 

0.035% 

1.265% 

0.002% 

1.706% 

0.002% 

0.321% 

0.401% 

0.993% 

0.010% 

0.808% 

0.007% 

0.878% 

0.237% 

0.000% 

0.000% 

0.000% 

(0.1,  0.3) 

0.015% 

1.607% 

0.012% 

1.670% 

0.683% 

3.348% 

1.754% 

0.000% 

0.096% 

1.760% 

0.007% 

0.329% 

0.000% 

0.000% 

0.000% 

0.000% 

0.000% 

0.701% 

0.000% 

0.000% 

0.301% 

2.992% 

0.059% 

0.000% 

0.031% 

1.228% 

0.006% 

0.917% 

0.244% 

1.332% 

0.443% 

0.199% 

0.069% 

1.691% 

0.079% 

0.703% 

0.100% 

0.716% 

0.876% 

0.869% 

0.084% 

5.300% 

0.055% 

3.276% 

0.003% 

0.404% 

1.844% 

0.000% 

(0.1,09) 

0.022% 

1.023% 

0.050% 

5.141% 

0.025% 

0.077% 

1.237% 

1.919% 

0.033% 

0.617% 

0.182% 

5.524% 

1.228% 

0.000% 

0.837% 

1.125% 

0.002% 

0.293% 

0.049% 

0.681% 

2.133% 

0.560% 

0.000% 

0.000% 

0.042% 

1.785% 

0.083% 

3.065% 

0.698% 

0.351% 

0.959% 

0.783% 

0.001% 

3.099% 

0.205% 

1.981% 

0.162% 

0.000% 

0.031% 

0.327% 

0.113% 

0.293% 

0.048% 

0.000% 

0.164% 

0.000% 

0.018% 

0.000% 

(0.1,  1.9) 

0.006% 

2.584% 

0.256% 

2.190% 

0.018% 

0.235% 

0.033% 

0.236% 

0.131% 

0.200% 

0.171% 

3.693% 

0.026% 

1.901% 

0.119% 

0.030% 

0.023% 

0.017% 

0.354% 

1.068% 

0.074% 

1.386% 

0.138% 

1.044% 

0.055% 

1.238% 

0.207% 

1.786% 

0.089% 

0.704% 

0.068% 

0.327% 
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(40,400) 

Variable 

Cost 

Machine  Cost 

( 10,  12) 

(10,15) 

( 10,20) 

( 10,40) 

(UB-LB)  / LB 

(UB-LB)  / LB 

(UB-LB)  / LB 

(UB-LB)  / LB 

1.382% 

1.454% 

2.019% 

1 .066% 

0.598% 

0.424% 

0.915% 

0.203% 

(0.1,  0.3) 

0.497% 

1.581% 

1.595% 

0.295% 

0.897% 

0.138% 

0.108% 

0.245% 

0.798% 

1.075% 

0.796% 

0.582% 

0.834% 

0.934% 

1.087% 

0.478% 

1.412% 

0.503% 

1.580% 

0.540% 

1 .369% 

0.131% 

0.596% 

2.008% 

0.60 

(0.1,  0.9) 

1.619% 

1.555% 

0.376% 

0.021% 

0.837% 

1.635% 

0.321% 

0.215% 

0.074% 

0.128% 

0.359% 

0.043% 

1.062% 

0.791% 

0.646% 

0.566% 

1.614% 

0.753% 

1.260% 

0.656% 

1.773% 

1.398% 

1.498% 

0.638% 

(0.1,  1.9) 

1 .398% 

1.024% 

0.725% 

0.305% 

1.143% 

2.074% 

0.151% 

0.198% 

1.397% 

0.675% 

0.489% 

0.367% 

1.465% 

1.185% 

0.824% 

0.433% 

1.080% 

1.063% 

1.333% 

0.645% 

0.015% 

0.408% 

0.869% 

0.256% 

(0.1,  0.3) 

0.115% 

1.006% 

0.018% 

0.267% 

0.036% 

0.430% 

1.665% 

0.062% 

0.157% 

1.272% 

0.128% 

0.264% 

0.281% 

0.836% 

0.803% 

0.299% 

0.656% 

0.552% 

0.917% 

0.781% 

6.128% 

3.457% 

2.630% 

0.219% 

3 

0.75 

(0.1,09) 

1 .667% 

1.425% 

2.976% 

0.391% 

1.522% 

1.253% 

1.872% 

0.782% 

0.944% 

0.126% 

0.374% 

0.764% 

2.183% 

1.363% 

1.754% 

0.587% 

2.657% 

0.815% 

0.567% 

1.256% 

1.520% 

3.035% 

0.621% 

0.993% 

(0.1,  1.9) 

1.415% 

4.781% 

0.828% 

0.513% 

1.191% 

0.776% 

1.642% 

0.215% 

2.842% 

1.403% 

0.635% 

0.707% 

1.925% 

2.162% 

0.859% 

0.737% 

3.222% 

2.583% 

2.498% 

1.449% 

3.397% 

2.630% 

1.203% 

1.242% 

(0.1,  0.3  ) 

2.001% 

2.135% 

0.779% 

0.253% 

3.224% 

2.222% 

1.272% 

0.517% 

2.981% 

4.616% 

2.131% 

0.666% 

2.965% 

2.837% 

1.577% 

0.825% 

1.375% 

1.906% 

1.414% 

0.246% 

1.987% 

6.750% 

0.884% 

0.524% 

0.90 

(0.1,  0.9) 

5.662% 

7.248% 

1.552% 

0.858% 

1.570% 

1.664% 

2.937% 

0.865% 

2.099% 

0.988% 

4.536% 

0.638% 

2.538% 

3.711% 

2.265% 

0.626% 

1 .309% 

0.574% 

0.162% 

0.015% 

1.153% 

0.998% 

0.443% 

0.737% 

(0.1,  1.9) 

0.946% 

1.362% 

3.033% 

0.383% 

0.545% 

0.971% 

0.626% 

0.345% 

6.427% 

1.255% 

2.125% 

0.554% 

2.076% 

1.032% 

1.278% 

0.407% 
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APPENDIX  B 

EXPERIMENTAL  STUDY  RESULTS  FOR  MULTI-PERIOD  CAPACITY 
ALLOCATION  SOLUTION  PROCEDURE 


The  results  of  the  experimental  analysis  discussed  in  Chapter  3 are  illustrated 
here.  The  numbers  on  the  upper  left  comer  of  each  table  indicate  number  of  periods, 
number  of  machine  groups,  and  number  of  operations,  respectively.  GAP  indicates  the 
gap  between  the  upper  bound  and  the  lower  bound  on  the  minimum  value  of  the  objective 
function.  CPU  time  indicates  the  processing  time  of  each  problem.  The  last  table  displays 
overall  means  of  ten  problem  structures. 
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Operating  Cost 

ON 

© 

CPU  Time 

0.0392  44.32 

0.0387  47.73 

0.0459  49.98 

0.0384  52.35 

0.0435  106.72 

0.0320  44.00 

0.0409  44.00 

0.0357  47.51 

0.0422  53.01 

0.0485  58.38 

0.0387  41.69 

0.0368  43.83 

0.0414  47.13 

0.0443  50.14 

0.0403  58.93 

0.0358  42.46 

0.0439  46.25 

0.0369  49.48 

0.0395  52.73 

0.0479  53.82 

GAP 

(0.1,15) 

CPU  Time 

0.0386  43.67 

0.0446  46.90 

0.0427  49.22 

0.0470  71.02 

0.0398  77.23 

0.0387  44.33 

0.0369  44.98 

0.0416  49.87 

0.0395  51.14 

0.0464  59.04 

0.0430  43.17 

0.0440  41.85 

0.0437  42.78 

0.0486  50.92 

0.0484  58.55 

0.0403  41.74 

0.0439  44.49 

0.0463  51.57 

0.0523  54.71 

0.0472  58.17 

GAP 

d 

CPU  Time 

0.0372  39.88 

0.0427  43.28 

0.0457  45.97 

0.0451  66.68 

0.0491  55.81 

0.0419  42.85 

0.0448  44.43 

0.0486  49.21 

0.0453  53.38 

0.0439  51.42 

0.0448  38.67 

0.0430  40.86 

0.0425  45.75 

0.0519  53.17 

0.0492  56.90 

0.0351  42.29 

0.0457  43.94 

0.0436  45.86 

0.0454  52.73 

0.0437  58.88 

GAP 

r- 

d 

d 

CPU  Time 

0.0336  37.62 

0.0485  39.88 

0.0445  46.52 

0.0476  47.79 

0.0469  52.18 

0.0348  39.22 

0.0405  42.96 

0.0496  44.55 

0.0502  50.04 

0.0517  53.66 

0.0383  43.44 

0.0486  44.49 

0.0521  45.86 

0.0408  54.32 

0.0473  54.10 

0.0325  44.43 

0.0474  42.84 

0.0473  45.26 

0.0403  40.65 

0.0377  52.18 

GAP 

(0.1  ,0.3) 

CPU  Time 

0.0428  38.17 

0.0498  39.54 

0.0462  45.86 

0.0466  45.21 

0.0493  57.45 

0.0524  45.92 

0.0434  43.66 

0.0572  48.94 

0.0511  46.25 

0.0496  56.30 

0.0540  35.75 

0.0466  41.52 

0.0452  48.17 

0.0437  47.78 

0.0528  54.76 

0.0726  32.85 

0.0504  38.94 

0.0459  48.99 

0.0524  50.91 

0.0545  52.01 
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